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^! Abstract 

By using the moment algebra of the Vlasov kinetic equation, we characterize the integrable 
Bloch-Iserles system on symmetric matrices [Bloch ct al. (2009)] as a geodesic flow on the Jacobi 
group Jac(]R^'^) = Sp(IR^") (DH(]R^'^). We analyze the corresponding Lie-Poisson structure by 
presenting a momentum map, which both untangles the bracket structure and produces particle- 
■ type solutions that are inherited from the Vlasov-like interpretation. Moreover, we show how 

^ \ the Vlasov moments associated to Bloch-Iserles dynamics correspond to particular subgroup in- 

clusions into a group central extension (first discovered by [Ismagilov, Losik & Michor (2006)]), 
which in turn underlies Vlasov kinetic theory. In the most general case of Bloch-Iserles dynam- 
ics, a generalization of the Jacobi group also emerges naturally. 
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1 Introduction 

Recent developments on kinetic equations of Vlasov type have shown how different integrable 
systems (both finite and infinite dimensional) may emerge from appropriate moment closures 
of the corresponding Vlasov dynamics [Holm and Tronci (2009b)]. For example, this is the 
case for two integrable PDE's in the Camassa-Holm hierarchy [Camassa and Holm (1993), 
Holm and Ivanov (2010), Holm and Ivanov (2011)]. The relation between these integrable sys- 
tems and the underlying Vlasov kinetic description resides in the fundamental fact that mo- 
ments are momentum maps which then preserve the geometric structure of the basic kinetic 
system [Holm and Tronci (2009b)]. For example, the Klimontovich particle solution is also 
a momentum map that is preserved in the case of the Camassa-Holm (CH) equation and 
this produces the particle description of the well-known peakon solution of the CH equation 
[Camassa and Holm (1993)]. Analogous arguments apply to a finite-dimensional integrable flow 
on the space of symmetric matrices, recently proposed by Bloch and Iserles [Bloch and Iscrlcs (2006)]. 

This paper uses the Vlasov approach to provide a full characterization of the Bloch-Iserles 
system in terms of its geometric Lie-symmetry properties. After reviewing the main properties of 
Vlasov-like equations, this Introduction continues by describing shortly how CH-type systems 
emerge by taking moments of the Vlasov distribution. Later, the attention is devoted to a 
special type of moment algebra initially appeared in the context of accelerator beam physics. 
This moment Lie algebra provides the fundamental understanding of the Bloch-Iserles system, 
whose main properties are also reviewed in the last part of this Introduction. 
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1.1 Kinetic equations of Vlasov type 

Vlasov-type equations govern the evolution of collisionless multi-particle systems that are far 
from their thermodynamic equilibrium. This is a typical situation in plasma physics, which 
is the natural context of the Vlasov equation for electrostatic (Poisson- Vlasov equation) or 
electromagnetic (Maxwell- Vlasov system) interactions. The dynamical variable is a probability 
distribution on phase space /(q, p, t), which retains the entire statistical information on the 
system. In the simple case of electrostatic interactions, the Poisson- Vlasov equation maybe 
written as 

relative to the Hamiltonian 

^(/) = ^y"/(q,P,i)(l|pf + t//(q,t))dqdp, with AUf = lfdp. (2) 

Here {•, -Ig^^ is the canonical Poisson bracket on M?'^ and Uj is the nonlinear collective potential. 

As for its Hamiltonian structure, equation (1) is well known to be a Lie- Poisson system 
on the Poisson algebra J^(M^'^) of phase-space scalar functions. More precisely, equation (1) is 
Hamiltonian (cf. [Gibbons (1981), Morrison (1982)]) relative to the Lie-Poisson bracket (see 
Appendix (A.l) and references therein) 

where the pairing (•,•) : x T{M.'^^) — )• M is the usual L^— pairing and the dual space 

J-*(M2fe) ^ Den(M2fc) 

is given by distributions on phase-space. 
The Lie group underlying Vlasov-type equations is often considered to be the group of 
symplectic transformations on phase space [Alarsden (1982)]. However, as pointed out in 
[Marden and Wcinstein (1981), Marsden ct al. (1983), Scovel and Weinstein (1994)], the rigor- 
ous version of this Lie group is given by the group of strict contact transformations 

Cont(M2'=+\0) = e Diff(]R2'=+i) I y;*0 = , 

where 6 = p ■ dq — ds is the contact one- form on M^fc+i and Difr(M2fc+i) ig the group of 
diffeomorphisms of M^'^+^. At the infinitesimal level, the identification between the Lie algebra 
cont(M2^+\ 61) = {X G X(M2fc+i) | £^0 = o} of strict contact vector fields (here £ denotes Lie 
derivation and X{M?^~^^) is the space of vector fields on M^'^"^^) and the Poisson algebra T{M.'^^) 
of phase-pace functions goes back to Van Hove's thesis [Van Hove (1951)], which contains a 
proof in coordinates. 

Remark 1.1 (The group of quantomorphisms) The above picture generalizes to the case 
of a symplectic manifold (M, fi) that is prequantizable. In this case there exists a circle bundle 
P — > M together with a connection 9 of curvature 0. Strict contact transformations on P (such 
that (p*6 = 9) are necessarily automorphisms of the principal bundle. This explains the common 
notation Cont(P, ^) = Aut{P,9) emphasizing the ^-preserving nature of these automorphisms, 
which are also well known as quantomorphisms. This is the usual setting in quantization theory. 
For Vlasov-type equations, the most interesting case is when the symplectic form is exact, e.g. 
M = T*Q. This picture is developed in the last part of the present paper. 
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Remark 1.2 (The Klimontovich solution) The Vlasov equation exhibits relevant geomet- 
ric structures that were recently reviewed in [Holm and Tronci (2009b)]. In particular, the 
Vlasov equation allows for a dual pair of momentum maps that is identical to the dual pair 
structure existing for Euler's fluid equation. The left leg of this dual pair is given by the well 
known Klimontovich solution 

/(q, p,t)= [ w{s) 5iq - Q(s, t)) S{p - P{s, t)) d^s , (4) 
Js 

where s is a point on a submanifold 5 C M^'^ of dimension m (with m < 2k), endowed with 
the volume form w{s)d"^s. The above expression is a momentum map which returns the well 
known Klimontovich solution (single-particle) [Klimontovich (1967)] when dim 5 = (see also 
[Weinstein (1983)]). 

1.2 Integrable Vlasov moment closures 

While the Poisson- Vlasov system mentioned in the previous section is of paramount importance 
in plasma physics, its underlying geometry lacks a geometric property which is common to many 
other Lie-Poisson systems, such as Euler's equations for fluids and rigid bodies. Indeed, the flow 
of the Poisson- Vlasov system is not geodesic since the Hamiltonian (2) contains a linear term 
in /. Then one is led to ask about the possibility of a purely quadratic Vlasov Hamiltonian. 
This question has been recently pursued in a series of papers [Gibbons, Holm & Tronci (2005), 
Gibbons, Holm & Tronci (2008a), Holm and Tronci (2009b)], which showed how several inte- 
grable systems emerge naturally by taking moments of a Vlasov distribution undergoing a 
geodesic flow on Cont(M^*^"*"^, 0). More particularly, a geodesic Vlasov equation reads as 

| + {/,e*/} = o, (5) 

where / = f{z,t) is the Vlasov distribution, while G * f denotes convolution with a kernel 
Q{z,z'). In this context, the momentum map properties of the moment operation have the 
feature of taking the geodesic flow on Cont(]R^'^"'"^, 0) to another geodesic flow whose properties 
depend on how moments are taken. 

The most celebrated moment method involves tensor fields of the type 

In the context of integrable systems, the above quantities first emerged in the study of the Ben- 
ney integrable hierarchy [Bcnney (1973), Gibbons (1981)]. These moments are called kinetic mo- 
ments and their momentum map properties were investigated in [Gibl:)ons, Holm & Tronci (2008b)] 
together with their Hamiltonian structure, which is Lie Poisson on the Lie algebra given by the 
Schouten bracket on symmetric tensor fields [Nijenhuis (1955)]. Kinetic moments were shown to 
provide a general understanding of integrable PDE's of Camassa-Holm type [Gibbons, Holm & Tronci (2008a)]. 
This is related to the fact that moments project the group Cont{T*Q, 6) underlying Vlasov-type 
systems to its subgroups, e.g. point transformations Diff((3) and their extended semidirect- 
product version Diff (Q) (s) J^((5) (here Q is the particle configuration manifold). For example, 
the (one component) Camassa-Holm equation is a geodesic equation on the diffeomorphism 
group Diff(M). These transformations are regarded as canonical transformations on the phase 
space = T*M that are generated by functions h{z) = h{q,p) that are (homogeneous) lin- 
ear in the canonical one- form pdq. This class of linear functions in the canonical momentum 
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is naturally associated to the first-order moment, whose dynamics can recover the Camassa- 
Holm equation as described in [Gibbons, Holm & Tronci (2005)]. This case correspond to 
= pp'e~l'?~'^ I/". An analogous argument holds for the two-component Camassa-Holm (CH2) 
equation [Gibl>ons, Holm & Tronci (2008a)], which is a geodesic flow on the semidirect-product 
group Diff (M) (s) J-'(M) . This Lie group is generated by phase space functions that are (inhomo- 
geneous) linear in the canonical one-form. These functions are naturally associated to both the 
zero-th and first order moments, whose dynamics returns the two component Camassa-Holm 
equation [CJhcn, Liu & Zhang (2006)] for an appropriate choice of Hamiltonian. Recovering the 
CH2 system from the geodesic Vlasov equation requires setting G = e~l''~^ ' + pp'e~^'^~'^ 

Besides the integrable CH and CH2 equations, the Vlasov approach has been also shown 
to apply to another finite-dimensional integrable system, which was introduced by Bloch and 
Iserles [Bloch and Iscrles (2006)] and will be described in the following sections. This system 
was shown [Holm and Tronci (2009b)] to arise from a quadratic Vlasov Hamiltonian upon taking 
another type of moments, which are different from the kinetic moments mentioned above. These 
alternative moments will be called statistical moments and are presented in the next section. 
The geodesic flow associated to the Bloch-Iserles system is completely understood only in one 
special case, in which this system describes a geodesic on the Lie group Sp(M^'^) of symplectic 
matrices [Bloch et al. (2009)]. However, more general cases are possible for the Bloch-Iserles 
dynamics and this paper aims to provide a complete characterization of their associated geodesic 
flows. 

Remark 1.3 (Momentum map solutions of moment equations) Another geometric prop- 
erty common to many moment equations is that they possess special solutions that are also mo- 
mentum maps. These solutions arise from the remarkable fact that the Klimontovich momentum 
map (4) survives the process of taking moments [Holm and Tronci (2009b)]. Then, for example, 
the Klimontovich solution yields the particle description of the peakon solution for both the CH 
equation and a modified version of the CH2 system [Holm, O Naraigii &: Tronci (2009)]. As we 
shall see, an analogous argument also holds for the Bloch-Iserles system. 

1.3 Statistical moments and their Lie-Poisson structure 

As mentioned in the previous section, the Vlasov moment interpretation applies to the Bloch- 
Iserles system [Holm and Tronci (2009b), Bloch et al. (2009)], provided moments are taken ap- 
propriately. In particular, this approach is based on the definition of statistical moments, whose 
corresponding Lie algebra was studied in detail by [Scovcl and Wcinstcin (1994)]. This sec- 
tion reviews the basic properties of these moments in preparation for their application to the 
Bloch-Iserles system, which in turn is briefly recalled in the next section. 

Upon denoting z = (q, p), the n— th statistical moment is defined by the following symmetric 
contravariant tensor on M^'^: 

1 r ™ 

^»n(/) := — / z®™/(z) dz G V "^ = 0,1,2,..., (6) 

where 'V' denotes the symmetrized tensor product (S)*^^™, and Xq is the (constant) total number 
of particles: usually one fixes the normalization Xq = 1, although we shall consider a generic 
value of Xq . 

Remark 1.4 (Statistical moments of accelerator beams) Statistical moments were intro- 
duced by [Channcll (1983)], in the context of beam dynamics for particle accelerators. In acceler- 
ator beam optics, the dynamics of the first three moments governs the linear dynamics of particle 
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beams moving through a magnet lattice. When nonhnear effects are also taken into account, 
higher order moments must be considered and the moment hierarchy must be truncated properly 
[Channell (1995)]. Moment (Casimir) invariants were classified in [Holm, Lysenko & Scovel (1990)]. 

We shall refer to the quantities (6) as simply 'moments', unless specified differently. These 
moments take the Vlasov Lie-Poisson structure (3) to the following Lie-Poisson bracket: 

oo / 

{F,G}{X)=: 

n,m=0 

where the moment Lie bracket is defined as 

and the pairing (•, •) is given by tensor contraction. Here, A^^'^ denotes the symmetric part of 
the tensor A, while JI is the Poisson matrix that can be chosen canonical for simplicity: 

Also, AB denotes the one-index contraction between covariant and contravariant tensors, that 
is {AB)Y-- = Aij,„kB^^^- (analogously for BA = {BA)^J^- = B^'^'-'Aiji,,). We may rephrase 
the above considerations in the following statement: 

Theorem 1.5 [Holm, Lysenko & Scovel (1990)] Let 

oo / n \ oo 

ST{V) :=0 V =:^STn{V) 

n=0 \i=0 J n=0 

be the Fock space of symmetric covariant tensors on the symplectic space {V,S), endowed with 
the Lie bracket (8). The map 

a : {ST{VU;-])^{HV), 5„^^l5„,Jz®" 
is a Lie algebra homomorphism. 

Remark 1.6 (Arbitrary constant symplectic structure) Although the remainder of this 
section uses the canonical symplectic form, so that J = — J^^, the following discussions in 
this paper often refer to an arbitrary symplectic structure J on M'^'^. Then, when J is not the 
canonical structure, the corresponding Poisson tensor is given in matrix form by 

N := -J"^ 

as it appears in the Lie bracket (8). The notation N is chosen to be consistent with [Bloch et al. (2009)]. 

Remark 1.7 (Conventions on factorials) In the expressions above one has to recall that 
[Sq, S'^] = [Sn, S'q] = Vn > 0, as it can be seen by expanding polynomials as 

^-^ n! ^-^ n] 

n>0 n>l 



dF dG 

dXn dXm. 



(7) 
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For the Lie algebra coadjoint operator, this means that ad^^j Xn = ad^^ Xn-2 = 0, as it is shown 
by relabehng the indexes in the moment bracket: 



{F,G} = Y,{Xn 



i+m— 2 



5F 5G 



5X„ 5Xr, 



(adxc'/xv Xk, 



n,k 



An these cancelations can be incorporated by simply defining (— n)! := 0, which is the convention 
adopted throughout this paper. 

From the discussion above one can easily recognize the graded structure of the moment 
algebra, with filtration 

[STn{V),SUV)] C STn+m-2{V) 

Then, one can see that the space 5*72(1^) dual to second order moments is a subalgebra, i.e. the 
algebra Sym(y) = V* y V* of quadratic forms on V , where we recall the notation V* V V* := 
(y* (g) V*Y^'^. The largest subalgebra is given by Sym(V)©y*©R and it wih play a central role 
in the remainder of this paper. In order to better understand the structure of the subalgebra 
Sym(y) © © R, one can restrict the moment Lie-Poisson bracket (7) to consider only the 
first three moments 

Xoif) = I /(z) dz, = I z/(z) dz, X2(/) = z^V(z) dz. 

Upon choosing the canonical symplectic form J = — J~^, this produces the Lie-Poisson bracket 

dF dG 



{F,G}, (X2,Xi,Xo)=Tr X2 



dX2 ' dX2 



1/ 



dF . dG dG BF , 
+ ^1 l^^l^ - l^^l^ (10) 



+ ^0 



dXi dX2 dXi 8X2 
dF f dG^'^ 



dXi \dXi 



For vanishing (Xo,Xi), the above Poisson bracket reduces to Lie-Poisson dynamics on the Lie 
algebra Sym(M^'^) of symmetric matrices on phase-space, with Lie bracket given by 

[5*2, 'S'2] jj = 5*2 JS'2 — »S'2 JIS'2 . 

This Lie algebra is generated by homogeneous quadratic polynomials on phase space, which are 

well known to arise as the Lie algebra of the symplectic group Sp(M^'^) [Guillcmin and Sternberg (1984)]. 

On the other hand, the space of inhomogeneous quadratic polynomials (i.e. non-vanishing 

(Xq, -'^i)) is known in quantization theory as the Lie algebra of the Jacobi group [Guillemin and Sternberg (1984)] 

Jac(M2'=) = Sp(R2*^) (D H(m2^'), 

where H(M^'^) is the Heisenberg group on M^'"' x M. The semidirect-product structure (de- 
noted by (D) of the Jacobi group has attracted a certain attention over the last decades in 
the quantization community. For applications of the Jacobi group in quantum mechanics see 
[Berccanu (2007a), Bcrccami ('2007b), Bcrccanu (2008), Bcrccanu (2009)]. Also, applications of 
the Jacobi group to signal processing were recently proposed in [Shuman (2003)]. 

Remarkably, the moment Poisson bracket (10) applies to the Bloch-Iserles system without 
substantial modifications, as shown in [Holm and Tronci (2009b)]. The present paper focuses on 
this particular integrable system thereby determining its underlying symmetry group in terms of 
the corresponding Vlasov moment formulation. This construction involves a quadratic moment 
Hamiltonian, whose associated geodesic flow will be characterized in the following sections. 



7 



Remark 1.8 (Phase-space averages and momentum maps) The moment method is used 
in kinetic theory to take averages of microscopic quantities. This process possesses a particularly 
interesting meaning in terms of momentum maps. Indeed, a general statement [Marsden and Ratiu (1999)] 
ensures that any equivariant momentum map J : M^'^ g* (generated by the symplectic action 
of a Lie group G with Lie algebra g) produces a momentum map : Den(]R^'^) — )> g* via the 
relation 

{J{f),0 = I /(q,p) (J(q,p),0 d'^qd*^?. 

Consequently, we can state that 

any matrix Lie-Poisson system can be realized as a Vlasov moment equation of some kind, 

as long as the G-action on M^*"' admits an equivariant momentum map. For instance, any G- 
action on M.^ can be lifted to a canonical G- action on M?^ = r*M'^', which in turn produces 
always an equivariant momentum map. A particular example of this class is given by the 
relation c7(/) = / q x p /(q, p) d^'qd^'p, which takes the Vlasov Lie-Poisson structure to the 
rigid body Lie-Poisson bracket on so* (3) (i.e. the dual space of antisymmetric matrices). All the 
above still holds upon replacing the vector phase space M."^^ by a finite-dimensional symplectic 
manifold {V, O) acted on by G, e.g. V = T*Q for some configuration manifold Q. 

It is now clear that statistical moments can be considered as arising from averaging of the 
momentum map 

J:V^ST{V), J(z) = jlz^4 

while kinetic moments arise in the same context from the momentum map [Gibbons, Holm & Tronci (2005), 
Gibbons, Holm & Tronci (2008a)] 

J:T*Q^ ST*{Q) ® Den(Q) , J(q, p) = {p®"}„gj, 5(q - x), 

where ST{V) are the contravariant tensors on the symplectic space V, while ST*{Q) denotes 
the space of covariant tensor fields on the configuration space Q. 

1.4 The Bloch-Iserles system 

The Bloch-Iserles system [Bloch and Iserlcs (2006), Bloch et al. (2009)] is the ordinary differ- 
ential equation on the space Sym(M"') of real n-dimensional symmetric matrices given by 

x{t) = [x{tf,n], (11) 

where N is a constant n-dimensional antisymmetric matrix and [•, •] denotes the matrix com- 
mutator. This system is known to be Lie-Poisson on the space Sym(]R") endowed with the Lie 
bracket 

[S, S']^ := SNS' - S'nS, for ah 5, 5' G Sym(]R"), (12) 

with Hamiltonian 

h{X) = \Tr{X'). (13) 

To see this, we identify Sym(]R") with its dual via the positive definite inner product {S, X) := 
Ti^SX), for S,X £ Sym(M"). Then, relative to this pairing, one obtains the infinitesimal 
coadjoint representation 

adc X = XSn - NSX. 
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Since 5h/5X = X, the (left) Lie-Poisson equation (61) associated to the Lie algebra (Sym(n), [ , ]n' 
and the Hamiltonian h recovers the Bloch-Iserles system (11). The system can also be viewed 
as an Euler-Poincare equation, associated to the Lagrangian 

(.{S) = \tt{S^). 

The Lie group underlying the geometry of the Bloch-Iserles system will be determined in 
Section 2 below. We now recall from [Bloch and Iserles (2006)] and [Bloch ct al. (2009)] the 
principal properties of the Bloch-Iserles system. 



Lie algebra structure. Suppose that the the antisymmetric matrix N has corank d. We 
can thus write 

" N 




N 



where N is a 2A; x 2A; invertible matrix, with n = 2k + d. In this case, upon denoting by A^dx2fc 
the space of d x 2/c matrices, we have the Lie algebra isomorphism 



^ : {{Sym{R^'')(i)Md>c2k] x^ Sym(]R'^), 



(Sym(R"), 



(14) 



given by 



with 



^{S,A,B) :-- 



S 

A 2B 



[(S, A, B), is', A', B')] c={[S, S'] ^ , AnS' - A'nS , (^N^'^) 



sym 



(15) 



see Proposition 2.4 in [Bloch ct al. (2009)] (note that we use here slightly different conventions). 

We now describe the Lie algebra structure of the left hand side of the isomorphism. The Lie 
bracket on Sym(M^''') is [•, •]^ defined in (12), with N replaced by N. The semidirect product 
Lie algebra SymeM."^^) M-dx2k is associated to the Lie algebra action of Sym(M^'^) on Mdx2k 
given by 

A^5-^:=-^N5, (16) 

where one verifies the condition [S, 5']^ ■ A = S ■ (S' ■ A) — S' ■ (S • A). This semidirect product 
Lie algebra is then extended according to the Lie algebra two cocycle 



C{iS,A),iS',A')) 



1 



AnA'^ - A'nA^) = {A^A'^) 



(17) 



where denotes the symmetric part of a matrix M. We denote by (Sym(]R^*'') (s)7V4rfx2fc) 

Sym(IR'^) the associated central extended Lie algebra and one can check that the cocycle (17) 
leads to the Lie bracket (15). Recall that a VF-valued two cocycle C on a Lie algebra q is defined 
as a bilinear antisymmetric map C : q x q ^ W such that 

C([^, V] , C) + cue, ^] , r?) + C([7?, C] , 6 = , V e, r/, C G 0. 

On the Lagrangian side, the isomorphism ^ yields the Euler-Poincare Lagrangian of the 
Bloch-Iserles system 



1{X) = ^ Tr{X^) = ^ Tr{S^) + Tr(^^^) + 2 Tt{B'^) =: £(5, A, B). 



(18) 



9 



We shall see later how this Lagrangian yields an invariant metric on an appropriate matrix Lie 
group. 

One can also find the corresponding Lie-Poisson Hamiltonian by simply considering the dual 
map to ^ in order to write the correct Lie-Poisson formulation on the dual Lie algebras. Upon 
using the trace pairing to identify Sym*(M'^) ~ Sym(]R"), we obtain the map 



given by 



^* : Sym*(]R") ^ Sym(]R2fc) x Mikxd x Sym(M'^) 

■{Y,K,M). 



Y K/2 



Indeed, we have the equalities 



((y, K, M), (5, A, B)) = Tr(y5) + Tv{KA) + Tr(MS) 



Y K/2 
K'^/2 M/2 

Y K/2 
K^ /2 M/2 



S A"^ 
A 2B 



,^{S,A,B] 



Thus, the Bloch-Iserles Hamiltonian (13) on Sym(M^'^) x M2kxd x Sym(M"') reads 
h{X) = ^ Tr(X2) = i Tr(y2) + ^ Tt{K^K) + ^ Tt{M^) =: h(Y, K, M). 



(19) 



We now consider two particular cases that are related to the moment dynamics associated to 
the Vlasov system. 



Case 1: N invertible. When n = 2k and the matrix N is non degenerate, the Lie bracket 
(15) restricts to its first component. In this particular case, the Lie algebra (Sym(2fe), [•, ■]^) 
can be identified with the Lie algebra of Hamiltonian matrices 

sP(m2*^, N"^) := {S e M2kx2k I + 5%"^ = O} . 

The latter is the Lie algebra of the symplectic group 

Sp(]R2^,N-^) = {5 G M2kx2k I /N"^^ = N"^} , 
and the Lie algebra isomorphism is given by 

r,p : (Sym(E2^-), [•, ■]^) ^ (sp(M2^N-l), [•, •]) , T{S) := N5 = S. (20) 

Thus, if N is invertible, the Bloch-Iserles system is the reduction of a geodesic flow of a left 
invariant metric on the linear symplectic group Sp(M2^N-^) [Bloch ct al. (2009)]. The Bloch- 
Iserles Lagrangian on 5p(M2fc,N-i) reads 

which corresponds to the following unreduced left invariant Lagrangian L : T Sp(M^'^, N~^) M: 
L{g,g) = Kg~'g) = ^Tr{^g{g)g) , with gig) = [N-'g-'f {N-'g-') = -g^N-^g. 
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In order to write the Lie-Poisson equations, we identify the Lie algebra sp(M",N~^) with its 
dual via the Killing form {X,S) = Tt(XS). One then obtains the coadjoint representation 
ad J X = [X,S] and the left Lie-Poisson equations 

[ ' 5X 

With respect to the same pairing, the dual map to the Lie algebra isomorphism F reads 

:sp(m2'=,N~^)* ^ Sym(M2*^)*, X ^ Vl^{X) = Xn = X. (21) 

Notice that, since the antisymmetric part of X'H vanishes, one has N^^A' + ^Y-^N"^ = 0, so 
that X G sp(M2*:,N-^)* ~ sp(]R2'=, N~^). Then, the Bloch-Iserles Hamiltonian (13) IS given on 
sp(M",N-i)* by 

h{X) = ^ Tr((<YN)^) = ^ Tr {X{-n^)X^) . 

Case 2: N of co-rank 1. A more interesting situation occurs in Bloch-lserles dynamics 
when d = 1, that is n = 2A; + 1 so that N is of rank 2k. In this case the Lie algebra isomorphism 
(14) becomes 

^: ((Sym(]R2'=)(s)(]R2fc)*^ XclR,[.,.]c) ^ (Sym(]R2'=+i), [., , (22) 

where the Lie bracket (15) is now written as 

[(52,5i,So),(5^,5;,5^)]^= ([52,5^]^, 5iN5^-5iN52, SiN^l^). (23) 

Notice that, for d = 1, the space of matrices Mdx2k reduces to row 2A;- vectors, rather than 
column vectors in M^'^. This is the reason why we have inserted the notation (M^^')* for the 
space of linear forms to distinguish them from ordinary vectors. Moreover, here we have chosen 
the suggestive notation {S,A,B) = (S'2,5i,S'o) to show that in the corank 1 case, we recover 
the Lie bracket underlying the dynamics of the first three moments of the Vlasov system, see 
(10). In terms of the variables 52, Si, So, the Bloch-Iserles Lagrangian (18) becomes 

e{S2, Si, So) = i tt{sI) + |5i|2 + 252. 

The expression (19) for the Bloch-Iserles Hamiltonian on Sym(M2^) x M^'^ x R becomes 

h(X2,Xi,Xo) = ^Tr(X|) + + ^Xl (24) 

Remark 1.9 (The geodesic Vlasov equation underlying Bloch-Iserles dynamics) Notice 
that the above Hamiltonian h{XQ, Xi, X2) can be written equivalently as a quadratic Vlasov 
Hamiltonian 

H{f) = l ||/(z)a(z,z')/(z')dzdz', with g{z,z') = ^ + ^z-z' + ^{z-z'f. 

Then, taking moments yields precisely the moment Hamiltonian h{XQ,Xi,X2). Thus, again 
this means that the Bloch-Iserles system possesses a Vlasov formulation in terms of equation 
(5). As we shall see, this allows to identify momentum map solutions of the BI system that 
arise from the Klimontovich map (4) . The problem of how the geometry of Vlasov dynamics is 
related to that of the BI system will be addressed in the second part of this paper. 
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Remark 1.10 (The role of the space Aidx2k) Consider the general case when N has corank 
d. One can use the isomorphism (20) (upon replacing N by N) to identify the Lie algebra 
(Sym(M2fc)©>j^^2fc) X Sym(M'^) with 

{5p(R^^,N~')®Mdy,2k) xcSym(R'^). 

If d = 1 we thus get the Lie algebra 

(sp(m2^',N~^)(s)(e2^')*) xcR. 

However, as we shall see in the next section, the case d = 1 suggests the need of incorporating 
ordinary vectors in M^'^ (rather than linear forms in (R^^)*) within the Lie algebra structure. 
We shall show how to slightly modify the isomorphism 

Sym(M2fc) (s) (R2fc)* ^ sp(M2'=, N"^) (S) (R^*^)*, (^2, Si) ^ (N^s, Si) 

into an isomorphism Sym(R^'^) (s) (R^'^)* sp(R^'^, N~^) ©R^*^, which naturally leads to the 
Lie algebra of the Jacobi group Jac(R2'', N~^) = Sp(R2'=, N~^) © H(R2'=, N~^). Consequently, in 
the general corank d case, we shall need an isomorphism that also affects the space Mdx2k- 

1.5 Summary of main results 

The results in this paper emerge from the discovery that the integrable Bloch-Iserles system 
possesses a Vlasov description in terms of statistical moments [Holm and Tronci (2009b)]. Then, 
this paper uses the polynomial algebra Po^M^*^) C J^(R^'^) associated to the Vlasov equation 
in order to identify the Jacobi group as the Lie symmetry underlying the integrable Bloch- 
Iserles (BI) system. We show that the general case when corank(N) > 1 requires a suitable 
generalization of the Jacobi group that is different from other generalizations already appeared 
in the literature. The Lie-Poisson structure of the BI system is then analyzed in detail and 
we show how this structure leads to a natural momentum map, which in turn arises from the 
Vlasov single-particle solution (Klimontovich) . 

The relation between Vlasov dynamics and the Bloch-Iserles system naturally leads to the 
question of how the Jacobi group is a subgroup of the strict contact transformations, which 
underlie Vlasov dynamics. This question is addressed in the second part of this paper, where 
explicit subgroup inclusions are presented. Also, we show how the generalized Jacobi group 
is a natural generalization of Cont(R^'^^^, 0). This particular result arises as a consequence of 
some recent developments on Abelian extensions of the quantomorphism group [Vizman (2010), 
Neeb and Vizman (2003)]. After the relation between the BI system and its underlying Vlasov 
description is completely characterized, the momentum maps accompanying the Klimontovich 
solution are analyzed in detail. 

2 The symmetry group of the Bloch-Iserles system 

As recalled above, when N is invertible the Lie algebra (Sym(R^'''), [•, ■]^) is isomorphic to the Lie 
algebra of the symplectic group Sp(R^'^, N~^). Therefore, the associated Bloch-Iserles system 
is the convective representation of the geodesies on Sp(R^'^, N~"^) relative to a left invariant 
Riemannian metric. More precisely, a curve g{t) G Sp(R^'',N~^) is a geodesic if and only if 
X{t) = g{t)~^g{t) solves the Bloch-Iserles system. 
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In this section we shall show that a similar geodesic interpretation exists for any antisym- 
metric matrix N. In order to do this, we will exhibit a Lie group whose Lie algebra is given by 
(Sym(M"'), [•, ■]^), where N is an arbitrary antisymmetric matrix of corank d. We first treat the 
case d = 1. 



2.1 Bloch-Iserles dynamics on the Jacobi group 

In order to clarify the occurrence of the symplectic structure J and its Poisson matrix operator 
N = — (see remark 1.6) as well as the role of the symplectic space M^''^ it is easier to use 
an intrinsic formulation rather than the matrix notation. We thus consider an arbitrary finite 
dimensional symplectic vector space {V, 0) with symplectic form Q. The corresponding notation 
is introduced in Appendix A. 2, to which the reader is addressed for further details. 



2.1.1 The Heisenberg and Jacobi groups 

The Heisenberg group associated to {V, il.) is the central extension of the symplectic vector space 
V (viewed as an Abelian group), relative to the group two cocycle ^ : V x V ^ M. given by 

3g{u,v) ■.= ^n{u,v), (25) 

so that the group multiplication is: 

(u, a){v, b) = + v,a + b + -^^(n, v) 

On the other hand, the Lie algebra is given by the central extension \]{V, J7) = F M of the 
vector space V, relative to the Lie algebra two cocycle ^. The Lie bracket on f)(y, M) is 

[(n, a), (?;,&)] = (0, (7(n, ?;)). 

Notice that the {left) Lie- Poisson bracket on the dual Lie algebra f)*(y, Q) = 1/* x M has the 
particular form 

{/, g}~ {k, m) = -mVL = -m {/, g}y, , (26) 

where { , }v* denotes the symplectic Poisson bracket induced on V*. This bracket is related to 
the Poisson bracket on V by the Poisson tensor : V* ^ V. For m € M fixed, the hyperplanes 
V* X {m} are Poisson submanifolds of f)*(y, fi). 

There is an action of the symplectic group Sp(y, Q) on H(y, Q) given by the representation 
on the first factor, that is 

{u,a) ^ g{u,a) := {gu,a), g£Sp{V,n). 

It is readily seen that Sp(y, 17) acts by group homomorphisms. Therefore, we can form the 
associated semidirect product 

Jac(y, Q) := Sp(y, 17) (S) R{V, n), 
called the Jacobi group. The group multiplication is given by 

{g,u,a){h,v,b) = (gh,u + gv,a + b+ ^n{u,gv) \ , 
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and the associated bracket on the Lie algebra }ac{V, Q) = spiV, Q) (S) f)(y, 0) reads 

[iS2,Si,So),iS'2,S[,S'o)] = ([52,5^] ,S2S[-S!,Si,niSi,S[)) . (27) 

In the special case V = M^, coadjoint motion on Jac(R^) has been studied in [Bcrndt (2003), 
Berndt (2006), Yang (2003)]. As we shall see in the next section, the Jacobi group (and its Lie 
algebra, above) fully characterizes the Bloch-Iserles system, in the case corankN = 1. 

2.1.2 Geodesic flows on the Jacobi group 

From the general setting above, we may now characterize the case corankN = 1 of Bloch- 
Iserles dynamics in a completely intrinsic fashion, that clarifies the occurrence of the symplectic 
structure J = — (see remark (1.6)) as well as the the role of the symplectic space M^'^ and 
its dual. Indeed, we recognize that the intrinsic formulation of the Lie algebra 

(Sym(M2'=)(D(]R2^)*) x^M 9 {S,A,B) 

is given by 

{Sym{V)®V*) xcR 3 {a,a,^), 

where 

Sym(y) = {a -.V xV a{u,v) = a{v,u), yu,v G V} 

is the space of symmetric bilinear forms on V. The Lie bracket is intrinsically written as 

[a, a']Yi = o" o o cr' — fj' o n" o fj, 

where we canonically identify an arbitrary symmetric bilinear form a : V x V ^ M with a 
symmetric linear map a'" : V ^ V* , via the usual relation 

(^a\u),v'j := a{u,v). 

For simplicity, we shall suppress the suffix b since the interpretation of a will be clear from the 
context. 

The semidirect product Sym(y ) © V* is taken with respect to the Lie algebra action 

a I— )■ (T • Q = —a o n'^ o 0" 
and the Lie algebra 2-cocycle C is given by the Poisson tensor: 

C{{a,a),{a',a')) =U{a,a'). 
Thus, the Lie bracket on (Sym(y) ©1/*) xc R reads 

[{a, a, /?), {a', a', /?')] = (^a a' - a' oU^ o a, aoU^ o a' - a a, U{a, a')) • (28) 

One can check that when V = M?'' and the matrix notation [11] = N is used, we recover the Lie 
algebra action (16), the cocycle (17), and the Lie bracket (23) given in the previous section for 
the case d = 1. 

The following lemma is a crucial result. It says that, for the corank one case, the Lie algebra 
underlying the Bloch-Iserles is isomorphic to the Lie algebra of the Jacobi group. 

Lemma 2.1 With the notation above, the map 

F: (Sym(y)(s)y*) Xc M ^ jac(y, J]), (a, a, ;3) ^ (^2, 5i, 5o) = (n« o a, n«(a), /3). 

is a Lie algebra isomorphism. 
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Proof. The linear map T is clearly an isomorphism. We now show that it respects the Lie 
brackets. Applying F to the Lie bracket (28) we get the expression 

(^n" o cj o o ct' - nf* o o nf* o a, U^{a o n» o a') - n^{a' o n" o a), n(a, a')) • (29) 

Using the formulas 

aoUKa = -a (u^{a)^ and Q (u^{a),U^{a')^ = U{a,a'), 

the expression (29) can be rewritten as 

(^n* oCTon^ocj' -n« otj'otf ofj, (tf o(j)(n»(a')) - (tf oa')(n«(a)), o(^n»(a),n«(a'))) . 

This recovers the Lie bracket (27) of the Jacobi Lie algebra, evaluated on the elements T{a, a, j3), 
T{a',a',P') oi]ac{V,n). ■ 

At this point, upon fixing 

{v,n) = (M2^-N^l), 

the Lie bracket on jac(M^*^, — N~^) yields the Lie algebra structure (23) for the Bloch-Iserles 
system (11) in the case corank(N) = 1. In this case, the above isomorphism is given by 

r : (^Sym(R2*^) (s) (R^'^)*) R ^ jac(R2*^, -N"^), (30) 

(52,5i,5o) ^ (52,5i,5o) = (N52,N5f,5o) . 
which produces the Bloch-Iserles dynamics on joc(R^'^, — N"^). 

Remark 2.2 (Nature of the Lie algebra isomorphism T) At this point, we observe how 
the isomorphism P : (Sym(y)(DF*) xc R jac(F, O), affects both the Sym(F) and the V* 
components in the Lie algebra (Sym(l/) ©y*) x^- R. Indeed, while the map 

Sym(y) 9 cj ^ n« o a € sp(y, O) 

was already introduced in (20), the map 

V* 3 a^nKa^V 

is needed here too. This fact is crucial for understanding the geometric construction underlying 
the Bloch-Iserles system and it is necessary in order to match the correct Lie algebra jac(y, Vt) = 

We can now rewrite the Bloch-Iserles Lagrangian on jac(R^'^, — N~^). We find 

£(52,5i,5o) = ^Tr(52^(-N-2)52) + 5f(-N-2)5i + 25o'. 

The dual isomorphism arises from the Killing form pairing {X,S) = Tt{XS) on sp(R^^, — N""*^) 
and it is given by 

joc(R2'=, -N"^)* ^ Sym(R2'=) x R^'^ x R 
(^2, Xu Xo) ^ {X2, Xi,Xo) = (^2N, -nx[, Xo) . (31) 
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This allows us to write the Bloch-Iserles Hamiltonian on the dual of the Lie algebra of the Jacobi 
group as 

h{X2, Xu^^o) = \ Tr(X|) + + ^Xl 

= i Tr {X2{-fe)Xl) + i^i(-N2)^f + ^-Xl 

where we used the expression of the Hamiltonian (24). 

Using the Euler-Poincare reduction theorem and the results obtained above, we obtain the 
following theorem. 

Theorem 2.3 Assume that the antisymmetric matrix N has corank(N) = 1. Then the Bloch- 
Iserles system X{t) = [X(t)^,N] describes geodesic motion on the Jacobi group Jac(2A:, — N~^) 
relative to the left invariant metric given at the identity by the Lagrangian 

£{So,Si,S2) = ^Tr(5j(-N-2)52) + 5f(-N-2)5i + 252. 

Proof. We know that, in the general case, the Bloch-Iserles system is an Euler-Poincare equa- 
tion on the Lie algebra Sym(M", [,]n)- When N has corank one, this Lie algebra is isomorphic 
to the Jacobi Lie algebra, by the preceding Lemma. The result follows by applying the Euler- 
Poincare reduction to the Jacobi group 3ac(M.'^'^, — N"^) relative to the left invariant Riemannian 
metric induced by the Bloch-Iserles Lagrangian. ■ 



2.2 Bloch-Iserles dynamics on the generalized Jacobi group 

Once we have characterized the special case corank(N) = 1), we are now ready to identify the 
Lie group symmetry of Bloch-Iserles dynamics for the most general case corank (N) > 1. This 
is the subject of the present section. 



2.2.1 Intrinsic definition of the Bloch-Iserles Lie algebra 



As we have seen, if N has corank d, then we split M" in such a way that the matrix representation 
of N has the form 

" N 




N 



where N is a 2A; x 2A; invertible matrix, with n = 2k + d. In general, if V is a vector space 
endowed with a antisymmetric (degenerate) bilinear form fi', then we can write V' = V (B W, 
where W = ker(r2') and V carries the symplectic form il. := il.'\v- As in (14), we have an 
isomorphism 

^ : ( (Sym(y ) (DL(y, W*)) xc Sym(I^), [■, ■]c) ^ (Sym(y W), [•, -Jn) 



given by 



^'(a,Q,/5) := 



a a 
a 2/3 



where we use the matrix notation to write a symmetric bilinear form on V' relative to the 
decomposition V' = V W. Here, L(y,W*) denotes linear maps V W* as usual. The 
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semidirect product Sym{V) (§) L{V, W*) involves the Lie algebra action of Sym(F) on L{V, W*) 
given by 

a a ■ a = -aoU'^ o a with a e Sym{V), a e L{V,W*) 
while the central extension by Sym(H^) is given by the cocycle 

C{{a,a),{a',a'))=U{a*{-),a'*{-)ry'- G Sym(VF), 

where n(Q;*(-), a'*(-))'^^™ denotes the symmetric part of the bilinear map 

{w,w')^U{a*{w),a'*{w')), w,w'eW. 

The Lie bracket on the left hand side is readily seen to be 

[{a, a, 13), {a', a', /3')] ^ = ([a, a']n, aoUK a' - a' oUK a , n(a*(.), q' *(-))^^"' ) , (32) 

where 

[a, a']n = o" o o cr' - cr' o n" o cj. 

The goal is now to obtain the generalization of Lemma 2.1 to the case corank(N) = d. We 
first define a generalization of the Heisenberg and Jacobi groups. 

2.2.2 The generalized Heisenberg group 

This generalization is given by the central extension 

R{V, Q; W) = L{W, V) x,^ Sym(W^) 3 {A, B), 
of the space of linear maps LiyV, V) from 14^ to ^ relative to the Sym(W^) -valued group 2-cocycle: 

^ : L{W, V) X L{W, V) ^ Sym(lF), ^{A, A!) := ^n{A{-) , A' i-)^'^ 

so that 

^{A,A'){wi,W2) = ^(^n{Aiwi),A\w2))+n{A{w2),A'{wi))y 
The group multiplication is thus given by 

{A, B){A\ B') = {A + A!,B + B' + SS{A, Al)) (33) 

and the Lie bracket is 

\{A,B)M.^')\ = (0,2^(A^O) = [^MA{-),Al{-W^) . 
If d = dim(VF) = 1 we recover the Heisenberg group since in this case L{W, V) = Y and 

m^AAl) = \^{A{^,A!{W^ = Ul{AA!). 

Notice that other generalizations of the Heisenberg group are also available; see e.g. [Yang (2002)]. 
However, to the authors' knowledge, the Lie group presented above appears in the literature for 
the first time. 
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Remark 2.4 (Lie-Poisson bracket on l)*{V,^l;W)) Upon choosing {V,^}) = (]R^^',JI) and 
W = R'^, we have L{W, V) = M.2kxd and the Lie bracket reads 

[{A,B),{M,B')] = (0, (^^M')'""") • 
We thus obtain the {left) Lie-Poisson bracket 

{/,5}-(/C,A1) = -Tr|^Al(^|^) (/C,>i) GA^rfx2fe xSym(M'^). 

This is the generahzation of the bracket (26) to the case where d>l. 

2.2.3 The generalized Jacobi group 

We now let the symplectic group Sp(y, 0) act on the generalized Heisenberg group H(y, 17; W) 
from the left by the action 

{AB)^{goA,B). 

Using the expression (33) for the multiplication in the general Heisenberg group and the property 
^^{gA,gA!) = ^{A,A'), for all g € Sp(U, 17), it follows that we have an action by homomor- 
phisms. Thus, we can define the group semidirect product 

Jac(U, Q; W) := Sp(U, n) (s) H(U, Q; W) , 

which we call the generalized Jacobi group. We now compute the associated Lie bracket. In 
general, for a semidirect product G(§)H oi two groups we have the formula 

In our case, the first factor is [S, S'] and the second factor is 

S • {A', B') - S' ■ {A, B) + (0, 2i^{A, A!)) = {S o A! - S' o A, n{A{-) , A! {■)yy'^) , 

we thus obtain the Lie bracket associated to the Lie algebra jac(y,r2;VF) of the generalized 
Jacobi group: 

[{S,A,B),{S',A!,B')] = {[S,S'IS o A! - S' o A,n{A{-),A!{-)yy'^) . (34) 

2.2.4 Geodesic flows on the generahzed Jacobi group 

We now formulate the generalization of Lemma 2.1 to the case corank(N) = d. 

Lemma 2.5 With the notation above, the map 

r : Sym(Uel^) ~ (Sym(U) (DL(y, VF*)) x Sym(VF) ^ jac(U, 17; VF), 

{a, a, /3) ^ (5, A, B) = {B^ oa,BK q*, /5). 
is a Lie algebra isomorphism. 
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Proof. We apply the imomorphism T to the Lie bracket (32). For the first component we 
obtain the desired expression B^oaoB'^ oa' — B^oa' oB^ oa = [S, S']. For the second component, 
we get 

BK{aoBK a'y - B^ o (a' o bK a)* = -B^ o a' o B^ o a* + BK a o bK a'* 

= SoA' -S' oA, 

which is exactly the second component of the Lie bracket on joc(y, fi; W). For the third com- 
ponent, we have 

S(a*(-),a'*(-)r"' = n{B^{a*{-)),B^{a'*i-))fy'^ = n{A{-), A' {■)f^'^ 

as required. ■ 

Then, upon fixing 

{V, n) = (M^'', -N"^) and W = R"^, 

we have 

H(]R2^\-N-1;M^) = M2kxd x,^Sym(]R'^) 3 {A,B) 

and the cocycle reads ^{A,A') = —^{A'^f^~^A'Y^^. Then, the Lie bracket on the generalized 
Jacobi Lie algebra jac(M^'^, — N~^;M'^) yields the Lie algebra structure (15) for the Bloch-Iserles 
system (11) in the case corank(N) = d. In this case, the above isomorphism is given by 

F: {Sjm{M.'"')®Mdx2k) Sym(M'^) ^ jac(M2^ -N"!; M'^), 

{s,A,B) ^ {s,A,B) = {ns,nA^,B) . 

Remark 2.6 (Specializations of the isomorphism F) Notice that, while the first compo- 
nent of F coincides with the isomorphism (20), another Lie algebra isomorphism 

r„ : Mdx2k xc Sym(M'^) ^ m^'', -N-^M'^) 
arises from the last two components of F. 

With the help of the isomorphism F, we can now rewrite the Bloch-Iserles Lagrangian on 
the Lie algebra joc(]R2*:, -N~^;R'^). We find 

£{S, AB) = ]^ TY(5^(-N-2)5) + TV {A^{-n~'^)A) + 2 Tr(B2). 

The dual isomorphism is given by 

jac(R2*^, -N~^ M'^)* ^ Sym(2A;) x Mafcxd x Sym(fi) (35) 

{X, /c, M) ^ (X, K, M) = {xn, -nic^M) , 

where {X, /C, M) € Sp(R2'=, -N"^) x Mdx2k x Sym(M'^). This allows us to write the Bloch-Iserles 
Hamiltonian on the dual of the Lie algebra of the Lie group Jac(2/c,d, — N~^) as 

h{X,IC,M) = lTT{X{-n^)X^) + ^Tr(/C(-N2)/C^) + ^Tr(A^2). (36) 
2 4 8 

Upon using the Euler-Poincare reduction theorem and the results obtained above, we obtain 
the following theorem. 
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Theorem 2.7 Assume that the antisymmetric matrix N has corank d. Then the Bloch-Iserles 
system X{t) = N] describes geodesic motion on the generalized Jacobi group 3ac(2k, d, —N^^) 

relative to the left invariant metric given at the identity by the Lagrangian 

This theorem concludes our discussion about how the Bloch-Iserles system describes a geodesic 
flow on the generalized Jacobi group. The next section will focus on the Lie-Poisson bracket on 
jac*(M^'^, — N~^;M'^) and will explore its underlying momentum maps. 



2.3 The Hamiltonian structure and its properties 

This section investigates the Lie-Poisson structure of the Bloch-Iserles system, starting from 
its underlying Jacobi group. We first compute the expression of the Lie-Poisson bracket on 
the Jacobi group and rewrite the Bloch-Iserles system in terms of the variables {X,IC,A4) G 
jac*(M'^'^, — N~^; M"^). Then we shall show how this Poisson bracket can be simplified (untangled). 
The latter process is made possible by the semidirect-product structure of the Jacobi group, 
so that reduction by stages [Marsden et al (2007)] (see Appendix A. 3) can be applied directly. 
The simplified Poisson structure is obtained by simple application of a Poisson diffeomorphism 
which takes the Poisson manifold jac*(M^'^, — N~^;M'^) into the direct-product Poisson manifold 
sp*(M^'^, — N^^) X f)*(]R^^, — N"-*^). Besides recovering an important class of Casimir functions, 
this process also provides us with a family of momentum map solutions, whose physical ex- 
planation will be clarified later in terms of the Klimontovich solution of the underlying Vlasov 
equation. 



2.3.1 The Lie-Poisson bracket 



We start our discussion by presenting the Lie-Poisson bracket associated to the generalized 
Jacobi group, which governs Bloch-Iserles dynamics. The (left) Lie-Poisson bracket on the dual 
Lie algebra joc(M^'^, — N~^;M'^)* is given by 



{/,g}_(Af,/C,M) 



{X,IC,M), 



6f 6f 6f 



6g 6g 5g 



Tr X 



5X' 5X 



-Tr LM 



5X' SIC' 6M J ' \6X' 6IC' 6M 
^6X6JC 6X6K: 



Tr JC 



(37) 



6)C J 




where {X,IC,M) £ sp{M?'', -N~^) x Mdx2k x Sym(]R'') and we used the expression (34) for the 
Lie bracket on the Lie algebra of the generalized Jacobi group. In few cases, it will be convenient 
to use the notation to distinguish between the left (gl) and right (g^) Lie-Poisson brackets 
accompanying the dual of an arbitrary Lie algebra g. 

Given a Hamiltonian h = h{X,IC,M), the associated Lie-Poisson equations are 



X 



' dh- 

^'Jx 



— o /C 



/C = /C 



6h 

Jx 



T 



(38) 



M = 
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where the diamond operator o : M2kxd x Mdx2k sp(M^*=, -N ^)* = sp(M^'^^ -N ^) is defined 
by the condition 

{AoJC,S) = {}C,SA) , for all 5 € sp(M2^ -N"^). 

We obtain 

{Ao}C,S) = {JC,SA) = Ti{JCSA) = Tr (^/CNN~^5) 

= Tr ((^/CNry'^N-^^) = ((^/CN)'^°'N-\5> , 

where we recalled that N~^5 is symmetric. 

If h is the Bloch-Iserles Hamiltonian (36) then evaluating the derivatives ^ = NA:'N and 
^ = — |N^/C-^ yields the explicit equations 

^ M = 0. 

which are the Bloch-Iserles equations on jac(M^'^, — N~-^,M'')*. 

As we have shown above, the Hamiltonian structure of the Bloch-Iserles system is the ordi- 
nary Lie-Poisson structure on the dual of the Lie algebra of the Jacobi group Jac(M2fc,_N-i;]R'^). 
However, using the special form of this Lie algebra in terms of its semidirect-product structure 

the next section will show that this Poisson structure can be simplified by using a Poisson 
diffeomorphism called the untangling map. This map can be naturally introduced in the context 
of Poisson reduction by stages, see [Marsden et al (2007)] and Appendix A. 3. This approach 
has the advantage of exhibiting the intermediate Poisson manifold r*Sp x f)* which is of the 
form T*G x P, where P is a Poisson manifold. In this case, there is a procedure to simplify 
(untangle) the reduced Poisson bracket. 

2.3.2 The untangling momentum map 

The Poisson bracket (37) possesses a special form, whose general setting is best described in 
[Krishnaprasad and Marsden (1987)] (see Appendix A. 4). As a consequence of this construc- 
tion, it is possible to transform the Lie-Poisson bracket (37) on jac*(M^'^, J; M"') into a direct-sum 
Poisson bracket on 5p*(M^'^, J) © fi*(M^'^, J; M'^). This property involves finding a particular mo- 
mentum map, which in turn produces interesting considerations on the system, such as Casimir 
invariants and momentum map solutions. As we shall see, the identification of this momentum 
map requires restricting f)*(M^^, J; M"') to an open subset. 

Prom the general theory (see [Krishnaprasad and Marsden (1987)] and Appendix A. 4), the 
main object is an untangling momentum map 

{)*(M2fc^j;M'i) ^sp*(M2fc,j) relative to the induced 
action of the symplectic group on f)*(M2'=,J;M'^) , which is given by 

iJC,M)^ilCg~\M). (40) 

From the abstract point of view, this is the representation of the Lie group G = Sp on the 
Poisson manifold f)*. This representation is induced by the homomorphism action of Sp on H, 
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which produces the semidirect-product structure Sp(S)H. The associated infinitesimal generator 
of 5 G5p(R2fc,j) is 

5^.(/C,7W) = (-/C5,0). 

Consider the Lie-Poisson structure on the dual of the generalized Heisenberg algebra f)(M^^, J; M'^) 

Since the symplectic group acts by homomorphisms, the action is canonical with respect to the 
Lie-Poisson structure. In order to show the existence of a momentum map we need to restrict 
the above bracket to the open subset det^A 7^ (cf. [Bloch et al. (2009)]). We shall denote 
this restriction by 

r(M2^JJ;M'^) = Mdx2k X s^iR'^) c JiM-^), 

where Sym(M'^) denotes the set of invertible symmetric d x d matrices (also known as 'inertia 
tensors'). We shall denote by the Lie-Poisson bracket restricted to the open subset 

^*(M2fc, j;]R'^). Also, we shall use the notation jac*(M2fc, j; M'i) = sp*(M2fc, j; M«f) x?i*(M2'=, J;R'^). 
A simple verification shows that 

Proposition 2.8 (Untangling momentum map) The map 

J^. :'^1(m2*^M;M'^) ^sp*(m2^J) 

(/C, M) ^ IC^M^^K. (41) 

is an equivariant momentum map relative to the action (40) of the symplectic group on the open 
subset ^*{M?^ of the generalized Heisenberg algebra. 

Proof. RecaU that the Poisson bracket on ^*(R2'^, J; M'') is 

{/.rfS.(C.>,):=-Iv(^M(|iJ,^). 

The momentum map property is verified by observing that 

{/, ( J^. , S) (/C, ) = - Tr (/C 5 = Tr (|^ (-/C 5)) = df ■ 5^,(/C, M), 
for all functions / = f{IC,A4), where the first equality follows from the formula 

The equivariance is verified as follows 

3^,{)Cg~\M) = ^r'g~^lC^M-'}Cg~' = ^gr'JC^M~'JCg~' = Ad;_i (j^,(/C,A1)) 
This completes the proof. ■ 
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Thus, upon writing the symplectic form J in terms of the Poisson matrix operator as 
J = — N~"^, we may untangle the restriction to jac (M^^, — N~"'^;M"') of the Lie-Poisson bracket 
(37), according to the following untangling map (see (64) in Appendix A. 4) 



u : {X,1C,M) ^ (X + Jr,{IC,M),IC,M 



X + ^NIC^M~^}C,}C,M ] =: {y,JC,M). 



Then, the untangled Poisson bracket { , on sp*(M^'^,JI) © f)*(M^'^, JJ; M*^) is the direct sum of 
the Lie-Poisson structures on sp*(M2'', -N~^) and 1)*{M.'^'', -N~^-i 



{f,9}uiX,IC,M) = -Ti\^X 
and the Hamilton's equations are written as 



SI 6g_ 
5X' 5X 



Tr M 



ST]' 

dJC J 



y 



y, 



6h 

sy 



1C = -M 
M = Q. 



6h 
6JC 



(42) 



See Appendix A. 5 for further details on the above equations. It is easy to recognize how the 
untangling map untwines the Lie-Poisson bracket, while it may entangle the Hamiltonian in a 
complicated way (see Appendix A. 5). However, this map is also important because it explicitly 
provides Casimir functions for the entangled bracket. Indeed, from the general theory (see 
corollary 2.3 in [Krishnaprasad and Marsdcn (1987)]) one has the following Casimir functions 
for the restriction of the Bloch-Iserles bracket (37) to jac*(M^''', -N^^;IR'^): 

1 



Cj{X,JC,M) = C.j{X + ^NJC^M-^IC 



where 



C,{X) = -TT{{Xf^] 



j = l...k 



is a Casimir function on sp*(M^'^,— N This result naturally recovers a class of Casimirs 
already found in [Bloch ct al. (2009)] by other methods. 



2.3.3 Momentum map solutions of the Bloch-Iserles system 

This section shows another consequence of the untangling momentum map in Proposition 2.8. 
Indeed, the equivariance property of this momentum map determines a special class of solutions 
of the Bloch-Iserles system. 

We first give the formulas for the adjoint and coadjoint action of the generalized Jacobi 
group. From the group multiplication law on Jac(M^'^, J; M*^), we obtain the following formula 
for the adjoint representation: 

Ad(3,^,6)(5',^',^') = (gS'g-\gA'-gS'g-'A,B' + {A^IgA'r^ - \ {A^IgS' g-^A^''^^ . 
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For the left coadjoint action we have 

Mlg j^t^y,{y,lC,M) = [gyg'^ [SAKg-^y^"^ - ]^AMA^S,lCg~^ -MA^S,M^ . 

Thus, there is a natural left action of the Jacobi group on the Heisenberg Lie algebra defined 

by 

{g,AB) ■ {A\B') := Ad(,,^,B)(0,^',S') = {0,gA',B' + {A^SgAT"') (43) 
whose dual left action reads 

{g, A, B) ■ (/C, M) = {JCg-^ - MA^S, M) . (44) 

These observations produce the following result: 

Proposition 2.9 The action (44) is canonical with respect to the Lie-Poisson bracket on the 
dual of the generalized Heisenberg Lie algebra. Moreover, when restricted to the open subset 
i)*[M?^,J;]R.'^), it admits the equivariant momentum map 

* :^;(m2'=,J;R'^) ^ jac*(]R2^,J;M'^) 

ilC,M) ^ Q r'^ IC'^M-^IC, /C, , (45) 

where the index + denotes the right Lie-Poisson bracket on f)*(M^'^,JJ;M'^). 



Proof. Since {0} x {}(M2fc, jj; M^) is a Lie subalgebra of jac(M2fc, j; M'^), the action (43) defines a 
Lie algebra isomorphism. Therefore, its dual action is canonical with respect to the Lie-Poisson 
structure on the dual. We now show the momentum map property. Since the infinitesimal 
generator reads 

(5, A, B)r^, (/C, M) = {-ICS - MA^I, 0), 

we get 

df ■ {S,A,B\^iJC,M) = - Tr (^^ICS^ - Tr (^^MA^?j . 
On the other hand, denoting ^p = (^', {S,A,B)), we have 

,/,(*,,.,A.„,.T..(>,(0,|^).-Tv(|^.M(|i)^,). 



where 

^ {^,{S,A,B))=Sr^lC^M~^ + A. 



die dJC 

The map ^'(/C, Al) = (^JI~^/C^A^~^/C, /C, A^) is therefore a momentum map. The equivariance 
follows from a direct computation. On one hand we have 

^{{g, A, B) ■ {KM)) = -^{Kg-^ - MA^^, M) 

= (\r^g-^1C^M-^lCg~^ - g~^ IC^ A^ ^ + \AKg-^ - \ AMA^^,]Cg-^ - MA^^,m\ . 
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On the other hand, we have 



Now we observe that these expressions are equal. ■ 

Remark 2.10 (Equivariance of ^) Note that one can add to ^ an arbitrary function of A4 
and stiU obtain a momentum map. In particular, the third component is an arbitrary function 
of M. However, in order to have equivariance, we need to make the above choice for ^. 

Being equivariant, the momentum map ^ is also a Poisson map, as long as both dual Lie 
algebras {)*(M^^, J; R"^) and jac*(M^'^, J; M'') are endowed with the (+) Lie-Poisson structure (or 
both have the (— ) Lie-Poisson structure). Therefore we have the following result: 

Corollary 2.11 Let h be a Hamiltonian on the dual Jacobi Lie algebra and define the collective 
Hamiltonian h := h o 'i'{JC',M') on 'i)*{R'^'',-'M''^;W^). Let K' and M' he solutions of the 
collective Hamilton's equations 




(46) 

forhon\)*{m?^,-n-^;W^). Then 

{X, JC, M) = (^-^ N'^ K'^M'-^K', /C', M'^ (47) 

is a solution of the Lie-Poisson equations (38) on the dual Jacobi Lie algebra. 

In the case of the Bloch-Iserles system, replacing (47) in the expression (36) we get the collective 
Hamiltonian 

h{lC',M') = ho^!{lC',M') = \Ti([K'^M'^^K'n^f)+-Ti{K'{-n^)K''^) + \Tv{M'^), 

8 V / 4 8 

and the Lie-Poisson equations (39) collectivize to 

t' = -- iK.' n^K'^M'-^ - M') /C'N 

2 ^ ' (48) 

, M' = Q, 
according to (46). 

The existence of momentum map solutions is a special feature of Vlasov moment equations. 
As explained in Remark 1.3, these solutions are inherited from the Vlasov dynamics through 
the corresponding Klimontovich solution. In order to understand the deep nature of this phe- 
nomenon, it is useful to study the various geometric structures that are shared between moments 
and their underlying Vlasov system. The next section addresses these questions for the special 
case of Bloch-Iserles dynamics on the Jacobi group and its generalized version. 
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3 Vlasov and the Jacobi group: subgroup inclusions 



Previous sections showed how the BI system is characterized by a well defined geometric con- 
struction based on the Jacobi group. As noticed in remark 1.8, the BI system also possesses 
a Vlasov formulation, which arises from a quadratic Vlasov Hamiltonian. In turn, the Vlasov 
equation possesses another interesting geometric footing in terms of Hamiltonian diffeomor- 
phisms and their central extensions (quantomorphisms, in particular). Then, it becomes a 
natural question to investigate how much of the geometry of the BI system is shared with 
Vlasov dynamics. As we shall see, the Jacobi group underlying BI is a subgroup of the quanto- 
morphism group underlying the Vlasov equation. Even in the most general case of BI dynamics 
(i.e. corankN > 1), we shall see that the generalized Jacobi group is a natural generalization of 
the quantomorphism group to the case of manifolds possessing a vector- valued symplectic form 
(see [Vizman (2010)]). More particularly, the present section shows that the moment closures 
determining BI dynamics correspond to particular subgroup inclusions into the group of quanto- 
morphisms of V xM. In the case when corankN = 1, these subgroups are given by M, Sp(y, $7), 
H(y, O), and Jac(y, 0). Their subgroup inclusions can be established by using an interesting 
group isomorphism between quantomorphisms and a central extension of the group of symplec- 
tic diffeomorphisms previously constructed in [Ismagilov, Losik & Michor (2006)]. The latter 
involves a group 2-cocycle, which we shall call the ILM cocycle. Moreover, the subgroup inclu- 
sions mentioned above can be extended to the general case of BI dynamics (i.e. corankN > 1), 
thereby showing that the generalized Jacobi group is itself a natural generalization of the ILM 
group cocycle. 

3.1 Background on prequantization central extensions 

Let (M, Q) be a connected symplectic manifold and suppose, for simplicity, that (M, M) = 
{0}. Denote by DifF(M, fi) = {ry € Diff(M) | r]*0, = 0} the group of symplectic diffeomorphisms 
of M (more precisely, we restrict to the connected component of its identity). This group is 
formally a Lie group whose Lie algebra X{M, 0) consists of Hamiltonian vector fields, that is, 
X(M, O) = {Xfi G X(M) I h € T{M)}. There is a natural central extension of the Lie algebra 
X{M, Q) associated to the exact sequence 

^ M ^ T{M) X(M, Q) 0, (49) 

where the Lie algebra brackets on T{M) and X(M, Q) are given by the Poisson bracket and 
minus the Jacobi-Lie bracket of vector fields, respectively. 

When O has an integral cohomology class, then there exists a principal circle bundle vr : 
P — )■ M over M and a connection one-form 6 € Q}[P) such that dO = 7r*J7. Such a symplectic 
manifold is said to be quantizable. In this case, the Lie algebra extension (49) can be integrated 
into a Lie group central extension 

^ Aut(P, 9) Diff(M, (50) 

of Diff(M,0), see [Kostant (1970), Souriau (1970), Ratiu and Schmid (1981), Vizman (1997)], 
called the prequantization central extension. Here Aut(P, ^) denotes the group {tp G Diff(P) | 
(p*0 = ip} of quantomorphisms (also called strict contact transformations) of (P, 6) (again, in 
what follows we restrict to the connected component of the identity) . Note that ^ is a contact 
form on P and that a quantomorphism is necessarily an automorphism of the principal bundle. 
In order to make clear the link between the central extensions (50) and (49), we recall that the 
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Lie algebra aut(P, ^) = {U G X{P) \ £1/6 = 0} of the quantomorphism group is isomorphic to 
J-{M), the Lie algebra isomorphism being given by 



U e aut(P, 6) ^ -iuO G T{M) 



(51) 



where iu9 G J~{M) is the function on M induced by the 5'^-invariant function iu6 on P. The 
inverse reads 



where C denotes the Reeb vector field, uniquely determined by the conditions i^^ = 1, i^^d^ = 0, 
and HorXf denotes the horizontal lift of the vector field Xf relative to the connection 6. For 
a trivial bundle P = M x M, we have Hor X = Xdm + a{X)ds and C = —ds- 

Remark 3.1 (Quantomorphisms vs automorphisms) When the condition = ^ is dropped, 
one is left with the group Aut(P) of automorphisms of P. In the general case, this larger group 
has no relation with Vlasov dynamics. However, it is interesting to notice that the Mawxell- 
Vlasov (and consequently Poisson- Vlasov) systems appeared in [Ccndra et al. (1998)] as coad- 
joint motion on Aut(M^'^''~^), at least for the Vlasov part of those systems. This was done upon 
introducing a redundancy in the construction of the system. In particular, this redundancy 
consists in adding the equations for particle trajectories, which in turn are already determined 
by the Vlasov equation alone. 

3.1.1 The ILM cocycle and its Lie algebra 

In the particular case when the symplectic form is exact, we can choose the trivial princi- 
pal bundle P = M x M. and the connection form 9 = a — ds, where = da. In this 
particular case, Aut(P, ^) is diffeomorphic to Diff (M,r2) x M, and the prequantization cen- 
tral extension can be described by a M- valued group 2-cocycle on Diff(M, fi), as shown in 
[Ismagilov, Losik &: Michor (2006)]. More precisely, given a point mo G M, the ILM cocycle, 
defined in [Ismagilov, Losik &; Michor (2006)], is 



where the integral is taken along a smooth curve connecting the point mo with the point 'r]2{mQ). 
The cohomology class of is independent of the choice of the point mo and the 1-form a 
such that da = 0, see Theorem 3.1 in [Ismagilov, Losik & Michor (2006)]. We denote by 
Diff (M,r2) x^,^^ M 3 {rj,a), the associated central extension of the group Diff(M, O), whose 
group multiplication reads 



Note that here we choose to work with the real line M instead of the circle S . The group 
isomorphism between the central extension and the group of quantomorphisms is given by 



/^HorX/-(/o7r)C, 



(52) 




(53) 



(r/i,ai) (772,02) = (r/i 0772,01 +02 + -Bmo(??i> %))• 



(77, a) G Diff (M, J7) xb M (^(^ G Aut(M x M, a - ds). 



(54) 



where ^(ri,a) is defined by 




(55) 



see [Ismagilov, Losik & Michor (2006)]. 
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Remark 3.2 (Point transformations) Among all Hamiltonian diffeomorphisms of an arbi- 
trary cotangent bundle M = T*Q, point transformations are regarded as cotangent lifts of diffeo- 
morphisms of the base configuration manifold Q. These transformations possess the fundamen- 
tal property of preserving the canonical one- form, i.e. a = p • dq on T*Q = R^'^. Therefore, the 
group cocycle (53) vanishes on the diffeomorphism group Diff((5), when the latter is interpreted 
as given by point transformations on T*Q. The same argument holds for the group of extended 
point transformations DiS{Q)(§)J-{Q). Notice that the groups T>iS{Q) and Diff (Q) (s) J^(Q) are 
the Lie groups underlying coadjoint orbits associated to the integrable Camassa-Holm equation 
and its two-component extension, respectively [Holm and Tronci (2009b)]. 

We now compute the Lie algebra 2-cocycle on X(Af , $7) induced by the ILM cocycle i?mo- 
Recall that this is given by the general formula 



C7^„(x,y)- I 



d 
ds 



s=0 



where r/j and sue curves in Diff(M, il) such that ^|^_Qr/t 



X and 4, 



ds \s=0 



6 



(56) 
Y. Also, 



recall that the Lie bracket on X{M, $7) xCm,, given by 

[iX,u),iY,v)] = i[X,Y],Cm,{X,Y)) . 
Lemma 3.3 The Lie algebra two cocycle associated to is given by 

Cmo{Xf,Xg) = {/q.^K^io) + {/,5a}("^o), 
where ha is the function on M defined by 

ha := h + a{Xh) ■ 

Note that we have dha = £xhOi. 

Proof. Let r]t and in Diff (M, 17) be curves tangent to Xj and Xg at t = 0. We compute 



d_ 

di 



d_ 

ds 



s=0 



B(,H.(.) = I 



i=0 



d 
ds 



s=0 Jrno 



6 (mo) 



(j/^a - a) = {£xfa,Xg) (mo) 



= {dfa,Xg) (mo) = {fa,g}{mo). 

Using formula (56) we obtain the required result. ■ 
The Lie algebra isomorphism associated to (54) reads 



aut(M X M, a - ds) ^ X(M, Q) x 



Cj7 



Xdm + kds ^ {X, k{mo)). 



Besides the Lie algebra isomorphism X{M,Q) xc^^ R ~ aut(M x M, a — ds), there is also a 
well-defined Lie algebra isomorphism with the space J~iM) of functions on M, 



h G T{M) ^ Xhdm + hads G aut(M x M, a - ds) 
obtained by particularizing (52) to our case. 



(57) 
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3.1.2 The case of symplectic vector spaces 

In this section, we specialize the above results to the case where M is the vector space V = M^'^ 
and is endowed with a constant symplectic form 0. We can suppose that we work in Darboux 
coordinates, that is, we have Q = dq Adp, z = (q, p) G M?''. We consider the central extension 
of Diff (y, for the particular choices 

mo := and a = ^i^' '^P ~ P ' ^^)' 

that is, we have 

rm{o) 

B{vi,V2) = / {r]la-a). 



Jo 

In this case, the cocycle C simplifies and we have the following 
Theorem 3.4 The Lie algebra two cocycle for {V, $7) and rriQ = is 

C{Xf,X,) = {f,gm. 
Proof. A direct check shows that for all / G J~iV), we have 

a(Xy)(z) = -id/(z).z and {aiXf),gm = -^{f,g}{0). 
Thus, applying Theorem 3.3, we have 

C{Xf,Xg) = {f + a{Xf),g}{0) + {f,g + a{Xg)}{0) 

= \{f,9m + l{f,9m = {f,9m- ■ 

The Lie algebra bracket on X(F, ft) xc M is thus given by 

[iXf,u),{Xg,v)] = {X{f^gy,{f,gm) , 

and the Lie algebra isomorphisms ^(V^) — >■ aut(y x M, q — ds) — )■ X{V,i}) x^ M simplifies to 
h ^ Xhdm + Kds ^ {Xh, /i(0)). 



3.2 Remarkable subgroups and Vlasov moment closures 

In this section we show how the groups M, Sp(y, O), H(y, 0), and Jac(y, 0) are naturally 
subgroups of the group of quantomorphisms of ^ x M, by exhibiting the concrete expression of 
group inclusions. For each cases, we compute the induced Lie algebra homomorphisms and we 
show that these subgroups correspond to moment closures of the Vlasov system on T{V)* . For 
each examples, we will use the Lie algebra isomorphism 

xcM^^(y), {Xf,u)^u + f - f{Q). 

As we shall see, the groups above are more naturally included in Diff (y, 0) x^ M rather than 
in the group Aut(y x M, ^) of quantomorphisms. 
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3.2.1 Translations on M 

There is a natural inclusion of the real line into Diff (y, 0) M given by 

R Diff(F,17) XijM, a^{id,a). 

This inclusion is a group homomorphism since (id, a){id, b) = {id, a + b + B{id, id)) = (id, a + b). 
The corresponding quantomorphism corresponds to translation by a on the M factor: 

(z, s) ^ f(id,a){^, s) = (z, s + a), 

since id = (see the notation ?] introduced in (55)). Taking the tangent map at the identity, 
we get the Lie algebra inclusion 

iR-.R^ X{V,n) XcM ^ T{V), u ^ {0,u) ^ {z ^ u), 

where the last expression means the function with constant value u, and the Lie bracket on M 
vanishes. The dual map to i^. is 

i*^:J'{V)*^R*, 4(/)= / /(z)dz, 

Jv 

and thus recover the moment Xo{f). 

3.2.2 Linear symplectic transformations 

There is a natural inclusion of the symplectic group into Diff (y, $7) x^ M given by 

Sp(F, 17) Diff (y, rj) Xb R, g^{g,0). 

This inclusion is a group homomorphism since, using -B|sp(y^f2) = 0, we have {gi,0){g2,0) = 
{gig2, B{gi, g2)) = {gig2,0). The corresponding quantomorphism reads 

(z,s) v3(g^o)(z,s) = (52, s), 
since ^ = 0. Taking the tangent map at the identity, we get the Lie algebra inclusion 

i,p : sp{v, n) ^ xiv, n)xcR^ Tiv), 

A ^ (z ^ Az, 0) ^ (^z^~^Tr {nAz®^)^ , 

where O denotes the matrix of the symplectic form, 5p{V,^}) is endowed with the usual Lie 
bracket of matrices, and J~^(y) is endowed with the Lie bracket given by the canonical Poisson 
structure. The last expression is obtained by noticing that z 1— )• Az is the Hamiltonian vector 
field associated to the function 

/(z) = Tr (l^Az®2) = i Tr (A^Uz^^) . 

Indeed, we have d/(z) = —z^QA and therefore Xf(z) = — (d/(z))''" = Az. The formula 

z,p(A)(z) = iTr(A^17z«2)^ 
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for the Lie algebra inclusion tells us that the appropriate isomorphism between Sym(V) and 
5p{V,Q,) is given by 5 = A^D., that is, 

r : 5 G Sym(y) ^ T{S) := -Q-^S = nS = Ae 5p{V, 0), 

where N is the Poisson tensor associated to il. We thus recover the isomorphism (20) used 
above. Identifying 5p{V, 0)* with 5p{V, Q) via the Killing form, the dual homomorphism is 

itp : :F{Vr^5piV,nr, i^^if) = -l^/(z)z^2^dz. 

Note that i*^ is consistently related to the second moment ^2(/) via the dual map F* : 
sp{V,n)* Sym(T/), r*(;f) = -Xn-^ = A'N, since 

r*oi*^{f) = x,{f). 

3.2.3 Heisenberg group 

Consider the natural inclusion 

R{V,n) ^mS{V,n) x^M, (v,a) ^ (Tv,a), 

where 

rv(z) := z + V 

is the translation. We will need the following lemma to show that this inclusion is a group 
homomorphism ([Ismagilov, Losik &: Michor (2006)]). 

Lemma 3.5 Consider the dijjeomorphisms t^^ and t^^ ofV. Then we have 

B{tvi,t^2) = ^^^(vi, V2). 

Using the lemma, we have 

(7"vi,ai)(rv2,a2) = (tv^ orvi,ai +02 + B{tv^,Tv^) = ^rvi+v2,ai + 02 + ^^^(vi,vi^ 

which proves that the inclusion is a group homomorphism. The corresponding quantomorphism 
reads 

(Z, S) ^ f{r^,a){2, S)= ^Z + V, S + a + ^^^(v, z)^ , 

where we used the identity 

P 1 
^v(z) = / «a -a) = B{t^, Tz) = -^^(v, z). 
Jo ^ 

We thus recover the natural action of the Heisenberg group on y x M. Taking the tangent map 
at the identity, we get the Lie algebra inclusion 

if, : f)(y, 0) ^ X{V, 0) xc K ^ aui{V x M, a - ds) ^ F{V), 
(w, n) I— 7> (z i-T- w, u) ^ w(?z + + — w^riz^ 5s I— > (z i-^- W"^r2z + ti) , 
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where we used that z i— )• w is the Hamiltonian vector field associated to the function /(z) = 
w-^Qz, and that we have the equahties /a(z) = /(z) — ^d/(z) • z = ^w'^Qz. Recah that the 
Lie bracket on l)(y,Q,) is 

[(v,n),(w,t;)] = (0,f](v,w)), 
and that on J-{V) it is given by the symplectic Poisson bracket. The dual homomorphism is 

i* : jr^vy ^ f,(y, ny, i;{f) = (-J^ /(z)z^i7dz, /(z)dz) 

The formula 

il^iw, u){z) = w'^fiz + u 

for the Lie algebra inclusion tells us that the appropriate identification between f)(y, Q) and the 
Lie algebra encoding the moments Xq , Xi is 

{Si, So)eV*xR^ (w, u) := j{Su So) = {-n~^Sf, So) € m 0) 

with dual map 

7* : {B,b) e f)(F,M)* ^ {Xi,Xo) = P*{B,b) = {n-^B^,b) eV xR. 
We consistently have the relation 

7*oi*(/) = (Xi(/),Xo(/)). 

3.2.4 Jacobi group 

Consider the inclusion 

Jac(T/, n) = Sp{V, n) (D H(1/, n) ^ Dm{V, Q)xbR, {g, V, a) ^ (tv o g, a). 
We will need the following lemma to show that this inclusion is a group homomorphism. 
Lemma 3.6 Given gi,g2 G Sp(y, il) and vi, V2 G V, we have 

B{t^^ ogi,T^2 °92) = ^ri(vi,5fiV2). 
Proof. Using the group 2-cocycle property of -B, we get 

-B(rvi 5l,Tv2 °g2) = B (Tvi ogi.Tva) +-B(rvi o gi orv2,g2) - B (Tv2,52) 

f>V2 /■91V2 /■V2 



/>V2 /■91V2 /■\ 

JO Jo Jo 

B (Tvi,TgiV2) - = -0(vi,5(iV2). 



Recall that the group multiplication on Jac(y, Q) reads 

(51, vi,ai)(52,V2,a2) = (gig2,^i + 5'iV2,ai + 02 + ^0(vi,yiV2; 
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We thus have 

(tvi o 51, oi)(rv2 o 02) = (tvi o 51 o ° 92, ai+a2 + B {t^^ o gi,T^^ o 52)) 

7"giV2+vi o 9192, ai +02 + ^0 (vi,5iV2; 

which corresponds to the group multipUcation on the Jacobi group. The injection is therefore 
a group homomorphism. 

The corresponding strict contact transformation is 

(Z, S) ^ (p(^^^og,a){z, S) = ^^Z + V, S + a + ^J7(v, ^z)^ , 

TV05(z) = / {g*T*a -a) = T*a = B{T^,gz) = -n{v,gz). 
Jo Jo ^ 

Taking the tangent map at the identity, we get the Lie algebra inclusion 

V, : }ac{V,n) ^ X{V,n) x^M ^ aui{V x M,a - ds) ^ J='(y), (58) 

{A, w, u) i-^' (z H> Az + w, u) {Az + w) 9z + (u + ^w^Jlz j dg 



smce 



^ ^z ^ ^ Tr (A^0z®2) + w^Oz + 



u 



where we recall that the Lie bracket on jac(y, 17) is 

[{A, V, u), {B, w, v)] = {[A, B],A-{w,v)-B- (v, u) + [(v, u), (w, v)]) 
= {[A,B],Aw - Sv,Jl(v,w)) . 

As before, the formula for the Lie algebra inclusion tells us that we need to consider the iso- 
morphism 

r:Sym(T/) x V* xR^}ac{V,n), 

given by 

{S2,SuSo) ^ {Aw,u) := {-n-'S2,-n-^sf,So). 

We thus recover the isomorphism (30). As we have seen, T endows Sym(y) x y* x R with the 
Lie bracket 

[{S2, Si, So), {S2,Si, So)] = {[S2, S2]m, SinS2 - S1NS2, Sinsf) , 

where N = —17^^ is the Poisson tensor associated to the symplectic form f]. The dual map 
r* : iac{V, fl)* — > Sym(y)* x 1/ x M consistently relates the dual homomorphism i*^^ : J^{V)* 
)ac{V,n)* 



C(/)(z) 



(^-^ ^ /(z)z«2^dz, - /(z)z^J7dz, /(z)dz 
to the three first moments {X2{f),Xi{f),Xo{f)), via the relation 

T*oil^{f) = {X2{f),Xi{f),Xo{f)). 
The main results of this section are summarized in the following theorem. 
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Theorem 3.7 The first three moments Xo(f), Xi(f), X2{f) of the Vlasov dynamics are as- 
sociated to the subgroup inclusion of the symmetry group of the Bloch-Iserles system into the 
symmetry group of the Vlasov dynamics, i.e. the inclusion of the Jacohi group into the quan- 
tomorphism group. These moments are given by the Poisson map obtained by taking the dual 
map to the associated Lie algebra inclusion. 

3.3 A commuting diagram underlying the Klimontovich map 

As we have just proven, there is a strict relation between the Vlasov kinetic equation and the 
Bloch-Iserles system when corank(N) = 1. We now show that, when corank(N) = 1, the momen- 
tum map solution of the Bloch-Iserles system obtained in §2.3.3 is naturally associated to the 
Klimontovich momentum map of the corresponding Vlasov equation. The latter is the solution 
provided by the single-particle trajectory and it is given explicitly by (cf. [Klimontovich (1967)]) 

f{z,t) = 6{-L-at)) =:iKiUC{t))- 

The map JKUm '■ V T{V)* is the momentum mapping associated to the natural left action of 
(r/, a) G Diff (y, fi) x^M on E 1/ given by C ??(C)- Being equivariant, JKHm is a Poisson map 
relative to the symplectic Poisson structure on V and the (+) Lie-Poisson structure on T{V)*. 
A review of the momentum map properties of the Klimontovich solutions was presented in 
[Holm and Tronci (2009b)]. Note that by applying the moment map / i-)- {X2{f), Xi{f), Xo{f)) 
to the Klimontovich solution we get the solution ^ (^'^ , ^, l) of the Bloch-Iserles system (in the 
case when Xo^f) = 1). This solution corresponds to a momentum map K2 : V — )• Sym(y)* xVx 
M*, whose last two components yield another momentum map Ki : V ^ V xM*. On the other 
hand, recall from §2.3.3 that there is another momentum map solution ^' : ^^(IR^fe, j; IR'^) ^ 
jac*(R^'^, J; M'^) of the Bloch-Iserles system, associated to the action (44) of the Jacobi group. 
We now show that, when d = 1, these momentum map solutions are naturally related by 
an equivariant momentum map Jh ■ V — > f)*(^i^^) associated to the natural action of the 
Heisenberg group on V. More precisely, upon recalling Remark 2.6, we will show that we have 
the following commutative diagram 

Sym(y)* X y X M* V xR* 



where 
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The upper left triangle commutes by the results obtained in the previous sections. As we just 
said JKiim is the momentum map associated to the left action of Diff (F, Q) x on V. When the 
subgroup Jac(y, 0) C Diff (F, Q,) x acts on V, by the subgroup action, i.e. {g, v,a)-C = sC+v, 
it produces the momentum map 

Jjac-.v ^jac*{v,n), 3jac{C) = (^-^CC^n,-c^n,iy 

Therefore, we have the relation Jjac = ijac ° ^KHm-, where i*^^ is the dual to the Lie algebra 
inclusion. Note that i*-^^ is itself a momentum mapping relative to the right coadjoint action of 
the Jacobi group on Tiy)* . Similarly, when the subgroup H(y, 0) acts on V by translations 
^ 1-^ ^ + V, it produces the momentum map 

and we have the relation Jj^ — Jftr/im? 

where ih is the Lie algebra inclusion. Note that both 
3h and 3jac are equivariant and hence Poisson with respect to the symplectic Poisson structure 
on V and the (+) Lie- Poisson structures. Besides the trivial relation 3h = i* o J jac between 3h 
and J Jac (here i denotes the Lie algebra inclusion i : f)(l/, f]) jac(y, fl)), we now show that 
J Jac = ^ o 3h- Indeed, we have 

^ o j^(c) = ^ (-c^j, 1) = Qj-i(-c^ji)^(-c^j), -c^j, 1^ 

This formula relates, in the case corankN = 1, the particular solutions given by the momentum 
map ^ and those obtained by the Klimontovich solution of the Vlasov equation via the first 
three moments. 

3.4 Remarks on the higher corank case 

While the previous sections addressed the question of how the Jacobi group is related to the 
group underlying Vlasov dynamics, it becomes natural to ask whether similar relations may ap- 
pear for the generalized Jacobi group, which underlies BI dynamics for corankN > 1. The first 
problem that we encounter in addressing this question is that there is no Vlasov moment opera- 
tion producing dynamical quantities in jac*(M^'^, J; M"^). Indeed, such moment quantities should 
arise from an averaging process involving a momentum map J : (M^'^,!) — t- jac*(M^'^, J; M"'), 
along the lines described in Remark 1.8. However, the generalized Jacobi group does not pos- 
sess a natural action on , J), so that the momentum map J is unlikely to exist. On the other 
hand, the generalized Jacobi group acts naturally on M2kxd, which in turn is not a symplectic 
vector space. Thus, we conclude that the case corank N > 1 of BI dynamics does not fit in the 
phase-space averaging principles outlined in the Remark 1.8. Consequently, one is led to seek a 
generalization of Vlasov dynamics which could include the BI system in the most general case. 
This section addresses the question of how one could develop a similar approach as above for 
the Bloch-Iserles system with an arbitrary corankN = d > 1. In this case, the resulting group of 
quantomorphisms (generalizing Vlasov dynamics) turns out to be the generalized Jacobi group 
itself. So, contrarily to the corank one case, there is no generalization of Vlasov dynamics that 
comprises the Bloch-Iserles system when d > 2. Or, in other words, this Vlasov dynamics 
coincides with the Bloch-Iserles system itself. 
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Recall that passing from d = 1 to an arbitrary corank d amounts in replacing the Ja- 
cobi group Jac(y, O) = Sp(y, fi) (S) H(y, O) by the generalized Jacobi group Jac(y, 0; W) = 
Sp(y, r2)(s) H(y,r2;VF) defined in §2.2.3. Contrarily to Jac(y, $7), the generalized Jacobi group 
is not a subgroup of Diff (V, M ~ Aut(y x M, a — ds). We thus seek a central extension of 

the Diff(y, generalizing the ILM central extension such that it has an automorphism inter- 
pretation and contains in a natural way the generalized Jacobi group. We already know that 
the cocycle involved should coincide with the cocycle of the generalized Heisenberg group, when 
restricted on Jac(y, Vl^). Therefore, the base space involved will be L{W,V), the extension 
is made by the vector space Sym(Ty), and the 2-form involved is no more a symplectic form: 
rather it is a closed vector-valued 2-form. Such an approach should thus be consistent with a 
generalized version of prequantization central extension by vector spaces. Such a theory indeed 
exists and has been recently developed in [Nceb and Vizman (2003), Vizman (2010)], as we now 
review very briefly. 

3.4.1 Abelian prequantization extension 

Let £^ be a vector space and uj € Q'^{M,E) be a closed i?- valued two-form on M. Suppose 
for simplicity that H^{M,E) = {0}, or, equivalently H'^{M) = {0} and let Diff(M,a;) be the 
connected component of the identity of the group of diffeomorphisms that preserves uj. When 
the period group of u is discrete, then there exists an Abelian principal bundle vr : P — )• M 
with structure group A = EjZ^^ and a connection one-form € il^(P, £■) such that = 7r*u}. 
In this case, we have a central extension of Diff(M, w) given by 

1 ^ ^ ^ Aut(P, 9) Diff(M, n) 1, 

as shown in [Necb and Vizman (2003)], where Aut{P,9) is the connected component of the 
identity of the group of automorphisms of the principal bundle that preserve the connection 
form 6. The associated exact sequence of Lie algebra reads 

0^ E-^ aut(P, 9) X{M, uj) 0. 

A very important point for the present discussion lies in the fact that the Lie algebra isomor- 
phism (51) is now replaced by the surjective linear map 

Ueaui{P,9)^ -kj9 G Tadm (M, E) , (59) 

where J-'admiM,E) denotes the space of admissible functions [Vizman (2010)] defined by 

J^ad,n{M,E) := {/ € T{M,E) \ df € Im(w^)} . 

Note that for all f,g& J-admiM, E), it is possible to define the bracket 

{{/, g}} := uj{X, Y) = df{Y) = -dg{X), (60) 

where X and Y are arbitrary Hamiltonian vector fields associated to / and g, respectively. This 
bracket is well-defined, since X and Y are determined up to a section of kei{u^). The second and 
third equalities are also well defined since d/lj^^^.j-^b-) = <ig\]<^cr{Lu^) ~ since / and g are admissible 
functions. It can be shown that when f,g€ Tadm'iM, E), then {{/, g}} G Tadm,{M, E) and {{ , }} 
is a Lie bracket on Fadm {M, E) . 

Let us observe that if the closed form u is exact, then the Abelian prequantization central 
extension can be defined by a global 2-cocycle, given by the same formula as the ILM cocycle 
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(53), although now taking values in E. We thus obtain a central extension Diff(M, cj) x E, 
with Lie algebra Xex{M,u}) ^Cmg ^ — ^{M,uj) xCmo ^' same formula (54) as above 

realizes a group isomorphism between Diff (M, cu) x b^^^ E and the group Aut(M x E,a — de) of 
quantomorphisms of the trivial bundle M x E, where de G i}^[E, E) is defined by de(n, v) = v. 



3.4.2 The space of admissible functions is finite dimensional 

Let us now consider the particular case when M = L{W^V), E = Sym(W^) and the closed 
Sym (VF)- valued 2-form on L{W, V) is defined by 

J7^(Ai?) := {A^nBY"'^. 

This is the natural generalization of the setting developed in §3.1.2 for the corank one case. We 
will however show that when d > 2, we get nothing else than the generalized Jacobi group itself. 

For an arbitrary (n x d) matrix A E L{W,V), we denote by Ai,i = its column. We 

thus have 

n^{A, B)ij = i {AfnB, + AjnBi) 

We denote by fij = fji : L{W,V) Sym(VK) the components of a function / : L{W,V) 
Sym(Ty), and by the partial derivative of / relative to the column Ak, so that dfij{A) -B = 

Ylt=i W^'^^' Using these notations we have the following result. (Notice that, in what follows, 
the summation convention on repeated indices is not used, unless otherwise specified). 

Lemma 3.8 The space of admissible functions associated to M = L(W, V) and the closed 
S'ym{W) -valued 2-form is 



TadmiL{W,V),Sym(W)) = / 



-—r- = 0, for all k ^ t,], — -p- = -— r^, tor all i ^ i > 
dAk dA^ 2dAj ^ ■' ] 

In particular, the kernel of Q}^ is trivial. 

- If d = dim(Vl^) = 1, there is no condition and J-admiY) = -^(^)- 

- If d = dim(14^) > 2, then the space of admissible function is finite dimensional, with dimension 
m(rn+i) ^ ^ _ jjj^^^^y^ _|_ (|iixL(M^) and consists of functions with matrix elements 

fij{Ai, Aj) = dj + ifi{Ai) + ipj{Aj) + B{Ai,Aj), 

where ipi, i = l,...,d are linear forms and B is a symmetric bilinear form. 

Proof. Recall that a function / is admissible if and only if there exists a vector field X G 
X{L{W, V)) such that n^{X{A),B) = df{A)-B, for ah A,B e L{W, V). Writing this condition 
for the matrix elements {i,j), we get 

lixrnB, + xJnB.) = ^^.B,. 

k=i ^ 

So, we immediately get the first condition = for all k ^ i,j, for all i,j. We also have 

= XfQ, for all i and = ^Xj'i}, for all i ^ j. This proves the second condition. 
It is clear that when d = 1 both conditions are empty and so all functions are admissible in 
this case, consistently with the statements in §3.1. 
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We now suppose d > 2. From the first condition, fij depends on the variables Ai and Aj 
only. Thus for i j, we get 

d'fij _ 1 ay,, _ d^n, _ 1 gy,, _ 

This means that /jj is necessarily of the form fij {Ai ,Aj) = Cij + ifij {Ai ) + tpij {Aj ) + Bij {Ai,Aj), 
where (fij and ipij are linear forms and Bij is a bilinear form. Using the fact that the second 
condition is valid for all i, we have ipij = (pj, ipij = ipi, and Bij = B. By symmetry, fij{Ai,Aj) = 
fji{Aj,Ai), which means (pi = ipi, Cij = Cji, and B is symmetric. We thus have 



dfji^^dfi^ 
dAj dA, 



2ipj + 2B{_,Aj 



and therefore fa are necessarily of the form fii{Ai) = Ca + 2ipi{Ai) + B{Ai,Ai). This proves 
the statement. ■ 



3.4.3 The generalized Jacobi group consists of generalized quantomorphisms 

Let us now assume that d > 2. Using Darboux coordinates on V, we can write elements in 
L{W, y) as Z = (Q, P), where Q, P are kxm matrices, 2k = n. We observe that all the results 
obtained in §3.1.1 regarding the ILM cocycle extend to the present case. It suffices to note that 
is exact, with potential form G il^{L{W, V), Sym(W)) given by 

a^(Q,P)(Q,P) = i(Q^P-P^Q)^'"\ 

Therefore, one can define the ILM cocycle exactly as before, and consider the central extension 
Diff {L{W, V),n^) X B Sym(T4^) with cocycle 

B{m,m)= / (??i*a^-a^). 
Jo 

This central extension is isomorphic to the group Aut {L{W, V) x Sym(H^), da^ — dS*) of au- 
tomorphisms of the bundle L{W,V) x Sym(T^) L{W,V) that preserve the Sym(l/K)-valued 
1-form — dS. This group plays the role of the quantomorphism group when d > 2. At the 
infinitesimal level, we have the Lie algebra isomorphisms 

X {L{W, V),n^) xcSym{W) ^ aut {L{W, V) x Sym(I^), - de) ^ TadMW, V), Sym(I^)). 

This shows, by Lemma 3.8, that the group Diff [L{W, V), x BS'ym{W) is finite dimensional. 
Note that the Lie bracket on the space of admissible functions is given by {{f,g}}, see (60). 

Of course, we still have the same group inclusions as in §3.2, for example Sym(VF), Sp{V, il), 
H(y, $7; W), and Jac{V, Jl; W) are included in the central extended group as, respectively, 

A^{id,A), g^{g,0), (V, S) ^ (rv, S), {g,Y,A)^{Tvog,A), 

where tv(Z) := Z + V is the translation. However, since the generalized Jacobi group and the 
ILM central extended group have the same dimension, we deduce that they are equal, being 
connected. Taking the tangent map to this Lie group isomorphism and composing with the Lie 
algebra isomorphism P defined in Lemma 2.5, we get the explicit identification 

Sym{V®W) ~ {Sym{V)(§) L{V,W*)) x Sym{W) ^ Tadm{L{W,V),Sym{W)), 

{s, A, B) ^ Q z^5z + {Azyy"" + B^= /(Z), 

between the Bloch-Iserles variables and the space of admissible functions. 
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4 Conclusions and open questions 



This paper showed that the Bloch-Iserles system on the space of symmetric matrices is a geodesic 
flow on the Jacobi group. This result arose essentially from the observation that the BI system 
possesses a Vlasov moment formulation, which in turn allows to identify momentum map solu- 
tions arising from the Klimontovich map. The Lie-Poisson structure of the general BI system was 
studied in detail and a special type of momentum map was used to transform BI on the Jacobi 
group Sp{V, Q) (D H(y, 0; W) into a system on the direct product group Sp(y, 0) x H(y, Jl, W). 

Moreover, this paper identified for the first time Vlasov-type systems with coadjoint orbits 
of the ILM group. Then, the Vlasov formulation of BI dynamics was shown to arise from 
a deep geometric construction, in which the Jacobi group underlying BI dynamics emerges 
as a subgroup of the quantomorphism (ILM) group. When trying to extend this result to 
the generalized Jacobi group, we found that the latter produces a natural extension of the 
quantomorphism group, so that the symplectic form may take values in the space of (invertible) 
symmetric matrices. This result was deeply inspired by the recent progress in [Vizman (2010)]. 

Many open questions emerge in the present work. Probably, the most intriguing concerns 
moments and its underlying Lie group. In [Scovcl and Wcinstcin (1994)], the authors mention 
(without giving any detail) that this Lie group is given by the jet group associated to quan- 
tomorphisms. However, this statement was never considered in later developments and this 
poses an important question concerning moment dynamics, which would then be explained as 
coadjoint dynamics on a group of jets. The same questions emerge for kinetic moments in §1.2. 
In this case, the Lie algebra properties associated to kinetic moments (dual to differential op- 
erators) [Gibbons, Holm & Tronci (2008b), Holm and Tronci (2009b)] are known to possess no 
underlying Lie group. Still, given the identification of kinetic moments with fiberwise polyno- 
mials, it is natural to wonder whether kinetic moments may allow for a geometric description 
in terms of a groupoid involving jets. However, this remains an outstanding open question that 
we plan to pursue in future work. 

A Appendix 

A.l Lie-Poisson and Euler-Poincare equations 

Recall that the ± Lie-Poisson bracket associated to a Lie algebra g is 



where ( , ) denotes the pairing between the Lie algebra and its dual. Given a Hamiltonian h on 
0*, the associated equations of motion are thus given by 



These equations are obtained by reduction of the canonical Hamilton's equations associated 
to a G- invariant Hamiltonian on the phase space T*G. In the formulas above, the upper sign 
corresponds to right invariance and the other to left invariance. When the Hamiltonian /i(/i) is 
regular (this is not the case for the Vlasov equation), the Legendre transform ^ = 5h/5fi takes 
to the Lagrangian side, thereby yielding the Euler-Poincare equations 





/i = =F ad^ fi. 



(61) 



d 5i 6e 



(62) 
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relative to a Lagrangian £ : g — t- M. These equations are also obtained by reduction of Hamil- 
ton's variational principle, associated to a G-invariant Lagrangian on TG. An Euler-Poincare 
approach to the Vlasov system has been proposed in [Cendra et al. (1998)]. Throughout this 
paper we shall consider both the Lie-Poisson and Euler-Poincare settings whenever the Legendre 
transform can be carried out. 

A. 2 Notation for symplectic vector spaces 

The symplectic form 17 : y x F — )• M generates the usual isomorphism 

n^:V^V\ {n\u),v) ■.= n{u,v) , 
whose inverse := is evidently given by 

n^:V*^V, ^{^\a),v) = {a,v) . 
In addition, the symplectic form generates the Poisson tensor 

n : y* X ^ M , n(a, /3) := Jl(!^«(a), ^l^{|3)) = {a, , 
which yields its associated Poisson bracket 

{f,g}{v)=U{df{v),dg{v)). 

Upon defining 

Ui:V*^V, n(a,/3) := (a, ntt (/?)), 

we obtain the relation 

The corresponding matrix relations among the symplectic form and its Poisson tensor may be 
easily found by fixing a basis on V, so that V* is endowed with the corresponding dual basis. 
We have 

where the square bracket notation indicates the associated matrix. Choosing a basis on V, we 
can write V = M."^^ and 

[n] = J, [n] = N 

with N = -J-^. 

A. 3 Semidirect-product reduction by stages 

Consider a semidirect product group G(g)H with group multiplication 

ig,h)ig',h') = igg',hig-h')), 

where g-h denote the action of G on H by group homomorphisms. This paper considers the case 
G = Sp(]R2fc, -N-i) and H = R{R'^'', -N'^-M.'^). The (left) Lie-Poisson structure on (gCDf))* is 
obtained by reduction of the canonical cotangent bundle T* (G ® H) by the whole group G(^H. 
This is the usual Lie-Poisson reduction. Note that this approach does not take into account of 
the particular semidirect product structure of the group. 
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However, one can also carry out the Poisson reduction in two stages, by making concrete use 
of the pecuhar form of the group. First we reduce the cotangent bundle by the normal subgroup 
H of G(S)i7 and we get the intermediate Poisson manifold 

T*{G(i)H)/H = T*G X f)* 

by the quotient map {g, g, h, h) i— )■ {g, g, h~^h) = {g, g, r]). The Poisson structure on T*G x t)* is 
the direct sum of the canonical Poisson structure on T*G and the Lie- Poisson structure on f)*. 
Then we reduce the Poisson manifold T*G x t}* hy G and obtain the dual semidirect product 
Lie algebra 

{T*Gx[)*)/G = (Q®l)r 

from the quotient map {g,g,'r]) i— ?> {g~^ ■ g,g~^ • r]). By the general result of Poisson reduction 
by stages ([Marsden et al (2007)]), the Poisson bracket obtained on (0(S)f))* by the reduction 
in two steps coincides with the Poisson bracket obtained by reduction in one step by the whole 
group G(g)H, that is, the Lie- Poisson bracket. 



A. 4 Untangling Poisson brackets 



Let (P, {•,-}p) be an arbitrary Poisson manifold and G a Lie group acting on the left on P 
by Poisson diffeomorphisms. Consider the Poisson manifold T*G x P endowed with the sum 
of the canonical Poisson bracket on T*G and the Poisson structure {•, -jp on P. Consider the 
diagonal action of G on r*G x P induced by the cotangent lift of left translation on G and 
the G action on P. As shown in [Krislinaprasad and Marsden (1987)], the reduced bracket on 
{T*G X P)/G Q* X P obtained by Poisson reduction is given by 



{f,g}ifJ;m) 



SI Sg_ 
6fj, ' 6fi 



+ {f^9}p + { dmg 



(m) 



dm /, 



Si 
5fi 



m)), (63) 



where the vector field on P denotes the Lie algebra action of ^ G g on P. We will refer to 
brackets of the form (63) as generalized Lie-Poisson structures. 

This type of Poisson structure arises in several physical contexts, such as the dynamics 
of charged fluids, superfluids and liquid crystals [Holm (1987), Holm and Kupcrshmidt(1982), 
Gay-Balmaz and Trouci (2010), Gay-Balmaz and Ratiu (2011)]. The remarkable fact about 
generalized Lie-Poisson structures is the existence of an untangling map, that is a Poisson 
isomorphism which takes the bracket (63) into a direct-product structure which does not explic- 
itly involve the g— action on P. Rather, this action is implicitly contained in the new dynamical 
variables. The only hypothesis for the existence of the untangling map is that the G-action on 
P allows for an equivariant momentum map [Krishnaprasad and Marsden (1987)]. Indeed, we 
recall the following general result 

Proposition A.l ([Krishnaprasad and Marsden (1987)]) Let G he a Lie group acting on 
a Poisson manifold P from the left. Suppose the G-action on P admits an equivariant momen- 
tum map Jp : P Q* and define the following Poisson structure on q* x P 

'SI sg; 

Sfj, ' Sfi 



{f^9}u 



+ {f,9}p, 



Then the map 



u:{q* X P, {■, •}) ^ (g* X P , {•, •}„) , (^, m) ^ {fi + Jp(m) , m) 
is a Poisson diffeomorphism. 



(64) 
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The name "untangling map" is due to [Holm (1987)], who showed how this map applies to 
electromagnetic charged fluids. In the present context, the special form of the Bloch-Iserles 
system allows us to untangle its bracket. Indeed, it suffices to consider the case 

T*G X P = T*Sp(M2^,J) X f)*(M2*=,J;M'^). 

As explained in Appendix A. 3 with the help of reduction by stages, the reduced bracket (63) in 
this case recovers the Lie-Poisson bracket on jac(M^'^, — N~^;M'^)*. 

A. 5 Alternative form of the Bloch-Iserles equations 

This Appendix shows the explicit form of the BI equations under the untangling map, that is 
under the change of variable 

One starts by writing the Bloch-Iserles Hamiltonian as 

h{y,K:,M) = ^Tr^(^3^ + ^N/C^7W-i/C^ (-N^) {y + M'^IC 

+ ^ Tr(/C (-N^) /C^) + ^ Tt{M'^) 
= ]^Ti{{ynf) - ^Tr(3^(-N2)/C^A^~^/CN) - ^ Tr (^(^(-N^) /C^A^-^/c)^^ 
+ ^Tr(^(-N2)/C^) + iTr(A^2). 
Upon using the formula 

one finds the explicit expression for the derivatives 

dy 2 

olC 2 2 

which, inserted in (42), produce the untangled Bloch-Iserles system in the form 

' y = [y,nyn] + ^ [y,'N'^JC^M~^icn] 
/C = -/C (NJ^N^)'^™ N"^ - ^ Jcn'^Jc^M-^m + ^ Mien 

[ M = 0. 
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